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ABSTRACT 


The  concept  of  a  limit  cycling  control  in  power  supplies 
and  speed  control  has  advanced  the  application  of  these  broad 
areas.   ON-OFF  switching  of  a  thyristor  provides  a  simple  and 
economical  method  of  control  and  regulation.   This  switching 
of  the  thyristor  causes  the  system  to  limit  cycle. 

Basic  analysis  and  design  of  speed  control  was  performed 
A  describing  function  was  developed  to  model  the  power-supply 
and  rectifier  bridge.   Then  it  was  used  to  predict  the  fre- 
quency and  amplitude  of  the  limit-cycle. 

A  digital  computer  program  was  used  to  simulate  the  sys- 
tem response  and  to  construct  the  describing  function  curves. 
To  verify  the  describing  function  validity,  the  limit  cycle 
predictions  were  compared  with  the  simulated  results. 

Fourier  analyses  were  performed  to  determine  the  ripple 
instability  and  subharmonic  effect  of  the  system  output. 
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I.   INTRODUCTION 

Limit  cycling  control  systems  have  enjoyed  considerable 
popularity  both  in  practice  and  in  the  literature.   The  prac- 
tical reasons  for  which  limit  cycling  systems  have  been 
developed  are  numerous,  not  the  least  important  of  which  is 
the  simplicity  of  an  on-off  power  element.   The  first  step 
toward  such  a  theory  was  actually  taken  by  Lozier  (Ref.  6)  at 
which  time  certain  qualitative  remarks  were  made.   Thereafter 
it  received  little  attention  In  the  English  literature.   The 
paper  by  Li  and  Vander  Velde  (Ref.  7)  represented  an  inde- 
pendent rediscovery  of  some  fundamental  principles  regarding 
limit  cycling  control  systems.   Work  of  a  similar  nature  has 
also  recently  been  published  in  the  U.S.S.R. 

In  this  discussion  stress  is  placed  on  considerations  of 
the  limit-cycling  system  as  a  control  system  and  its  signal 
transmission  properties  are  derived  as  well  as  the  loop  limit 
cycle . 

The  discussion  is  limited  to  negative  feedback  configu- 
rations with  a  nonlinearity  in  the  forward  loop,  such  as 
shown  in  Figure  1.1.   Here  N  stands  for  the  nonlinearity,  and 
G(s)  for  the  loop  linear  elements.   It  is  assumed  that  this 
is  a  limit-cycling  system.   That  is,  in  the  absence  of  an  in- 
put, the  system  sustains  a  nearly  sinusoidal  undamped 
oscillation  at  its  output.   If,  as  is  commonly  the  case,  the 
function  of  this  system  is  to  reproduce  the  input  signal  at 
a  higher  power  level  at  its  output,  then  its  response  to  a 
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FIGURE  1.1:   Simple  limit  cycling  control  system. 


r(t)  input 


T^V     Ave-RAGE  Output 


Actual    0utput    C  (0 


FIGURE   1.2:      Typical'  input    and   its    corresponding   output 


FIGURE   1.3:      Model   of   x(t)    over   any    limit    cycle. 


slow  input  signal  might  look  as  shown  in  Figure  1.2.   For  a 
well-designed  control  system  the  output  follows  the  input  on 
the  average,  within  some  dynamic  following  error.   Over  any 
period  of  the  limit  cycle,  T  ,  an  approximate  model  of  this 
error  signal  is  shown  in  Figure  1.3.   In  this  model  the  sin- 
usoid of  amplitude  B  is  associated  with  the  limit  cycle,  and 
the  D-C  bias  of  amplitude  a  with  the  input  forced  following 
error . 

In  this  study  a  system  which  has  properties  mentioned 
above  will  be  examined.   The  specific  system  to  be  examined 
is  an  armature  voltage  controlled  D-C  motor.   Its  armature  is 
supplied  with  a  poly-phase,  thyristor,  A-C  to  D-C  rectifier 
and  its  field  excited  by  a  separate  D-C  source. 

The  goal  in  this  study  is  to  investigate  the  control  of 
speed  of  this  D-C  motor,  operating  under  forced  limit  cycling 
conditions.   A  mathematical  model  is  proposed  and  then  ana- 
lyzed using  frequency  response  techniques. 

The  thyristor,  which  is  used  in  the  rectifier  bridge  is 
a  four  layer,  three  terminal  solid-state  device.   Functional- 
ly it  is  equivalent  to  its  gas  tube  counterpart,  a  thyratron. 
The  characteristics  and  circuit  symbol  of  the  thyristor  are 
shown  in  Figure  1.4.   In  the  forward  direction  the  device  has 
two  stable  modes  of  operation,  one  of  high  impedance,  of  sim- 
ilar amount  to  that  in  the  reverse  direction  and  the  other, 
of  very  low  impedance.   In  general  the  forward  break-over 
voltage  is  inversely  proportional  to  the  gate  current.   This 
indicates  that  the  device  could  be  controlled  by  controlling 
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FIGURE  1.4:   Thyristor  transfer  characteristics 
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the  level  of  gate  current.   The  control  of  the  forward  break- 
over voltage  provides  the  interval  of  the  forward  conduction 
and  thus  the  amount  of  power  per  cycle  may  be  regulated.   In 
this  paper  it  is  assumed  that  the  entire  rectifier  bridge  is 
controlled  as  a  unit  rather  than  by  individual  device  for 
each  thyristor.   The  gates  of  all  thyristors  in  the  poly- 
phase bridge  are  connected  in  parallel  and  a  single  firing 
signal  is  provided  to  all  devices.   A  simplified  circuit 
diagram  is  contained  in  Figure  1.5. 

Conceptually  the  system  operates  as  follows:   assume 
that  the  output  speed  is  below  the  reference  speed.   The  rec- 
tifier bridge  is  then  supplied  with  an  ON  signal,  and  the 
bridge  output  to  the  armature  of  the  D-C  motor  is  a  three- 
phase,  full-wave  rectified  voltage.   The  filter  action  of  the 
motor  alternates  the  harmonics  of  this  voltage  and  provides  a 
relatively  smooth  speed  at  the  output.   When  the  speed  rises 
above  the  reference  level  the  bridge  is  supplied  with  an  OFF 
signal.   There  is  a  time  delay  before  the  bridge  turns  off, 
because  the  shut  off  conditions  for  all  thyristors  must  be 
met.   The  time  delay  is  approximately  of  the  order  of  h 
cycle  of  the  supply  frequency.   In  the  OFF  state  the  input 
voltage  to  the  armature  of  the  D-C  motor  is  zero  with  reverse 
current  essentially  blocked.   The  bridge  will  remain  in  the 
OFF  state  until  the  output  speed  falls  below  the  reference, 
whereupon  the  cycle  will  repeat  itself.   This  limit  cycling 
of  the  system  is  the  particular  area  to  be  considered  here. 
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FIGURE   1.5:      Simplified   circuit   diagram   of   the    system 
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II.   BASIC  ANALYSIS  AND  DESIGN  OF  SPEED  CONTROL 

D-C  motors  with  shunt  excitation  are  particularly  appli- 
cable for  speed  control  systems.  The  speed  of  this  motor  can 
be  expressed  by 

n  _  Va-iaR 

^  "  -CT$ (2_1) 

where, 

a)   =  Speed  of  motor 

Va  =  Armature  voltage 

ia  =  Armature  current 
R  =  Resistance  of  armature  circuit 
$  =  Operating  magnetic  flux 

Ci  =  Coefficient  of  proportionality 

It  is  apparent  from  equation  2-1,  that  the  speed  of  such 
motors  can  be  controlled  either  by  varying  the  armature  volt- 
age, Va,  or  the  operating  flux,  $,  that  is,  by  varying  the 
field-excitation  current  of  the  motor.   However,  the  advan- 
tage of  speed  control  characteristics  of  a  D-C  motor  cannot 
be  fully  utilized  unless  an  adjustable  D-C  voltage  can  be  ap- 
plied to  the  armature.   In  order  to  take  full  advantage  of  a 
D-C  motor  characteristics,  it  is  necessary  to  provide  an  ad- 
justable output  voltage  converter  which  would  convert  the  A-C 
line  voltage  into  an  adjustable  D-C  voltage,  which  can  be  ap- 
plied to  the  motor. 
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The  very  popular  rotating  converter  consisting  of  an  A-C 
motor  driving  a  D-C  generator,  and  known  as  the  "WARD  LEONARD" 
system,  is  one  of  the  most  widely  used  adjustable  voltage 
supplies. 

The  rectifier  converter,  combined  with  complete  elec- 
tronic control  and  regulating  system  is  perhaps  the  most 
practical  unit  available.   Its  prime  advantage  over  the  ro- 
tating type  converters  is  that  only  one  rotating  machine  is 
used  instead  of  three.   In  addition,  the  compactness  of  the 
electronic  system,  the  flexibility  for  adaptation  to  various 
conditions,  accuracy,  and  sensitivity  are  improved. 

The  particular  interest  here  is  to  control  the  speed  of 
a  D-C  motor  by  controlling  the  conduction  angle  of  the  thyris- 
tors,  which  are  used  as  an  A-C  to  D-C  converter  to  supply  the 
armature  of  the  motor  as  shown  in  Figure  1.5.   In  this  config- 
uration a  feedback  control  system  was  employed  to  achieve  the 
desired  performance. 

The  system  in  its  simplest  block-diagram  form  is  illus- 
trated in  Figure  2.1.   In  order  to  permit  analysis,  the  system 
components  are  lumped  and  partitioned  as  illustrated  in 
Figure  2.2.   In  this  figure,  block  "G2"  represents  the  motor 
and  load.   It  is  assumed  that  they  are  entirely  linear. 
Block  "G-l"  is  included  to  take  into  account  any  gain  that  may 
be  associated  with  the  comparator,  firing  circuit,  and  com- 
pensation network  which  may  be  added  to  the  system.   Again, 
this  block  is  assumed  to  be  linear.   The  source  and  rectifier 
bridge  are  represented  as  block  "N"  which  is  assumed  to  have 
a  periodic  input  and  output  for  development  of  basic  describing 
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function  (which  will  be  mentioned  later  in  this  study).   The 
general  nature  of  this  block  is  nonlinear.   Block  "H"  repre- 
sents a  conversion  factor  of  the  speed  to  the  voltage.   It 
can  be  considered  as  a  feedback  gain. 

The  principle  of  operation  can  be  described  as  follows: 
the  speed  setting  voltage,  Vr,  represents  the  desired  motor 
speed,  and  the  tachometer  gives  K^  per  radian  per  second, 
which  is  directly  proportional  to  the  tacho  generator  speed. 
Regulation  of  tachometer  voltage  is  equivalent  to  the  regula- 
tion of  motor  speed.   The  difference  between  Vr  and  the 
tachometer  voltage,  Ktoom,  is  an  error  signal  which  can  be  rep- 
resented as 

a  =  Vr  -  Ktum  (2-2) 

This  signal  shows  the  deviation  of  the  motor  speed  from  its 
desired  value.  The  purpose  of  the  automatic  speed  control  is 
to  minimize  this  error  signal  and  obtain  the  desired  value  of 
the  speed.  This  error  signal  operates  the  firing  circuit,  it 
controls  the  motor  speed  by  adjusting  the  conduction  angle  of 
the  thyristors  in  the  rectifier  bridge,  thus  altering  the 
average  current  and  voltage  value  of  the  armature  circuit. 

One  procedure  for  analysis  of  the  system  is  by  using  the 
conventional  transfer  function  concept.   The  system  transfer 
function  is 

fi(s)  _     Gi (s)G2(s)N 


Vr(s)    l+Gi(s)G2(s)H(s)N  (2-3) 

The  linear  portion  of  the  system  may  be  block  diagrammed 
as  shown  in  Figure  2.3- 
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FIGURE  2.3:   Block  diagram  for  linear  system 
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The  transfer  function  for  the  motor  and  load  illustrated 
in  Figure  2.6  is  developed  as  follows.  The  loop  equation  for 
the  armature  circuit  is: 

Va  =  iaR  +  L  ^+  KmWm         (2-4) 
Torgue-equilibrium  equations  are 

Td  =  Kt  ia  (2-5) 

Ti  =  Jwm  +  f  a)m  (2-6) 

From  equations  2-5  and  2-6 


0)m   =  -_^_  o)m  +  _±_  ia 


and  from  equation  2-4 


R  ,     K   ,.  ,  1 
L  ia  -  L"  wm+  r 


ia  = r  xa  "  f—  wm+  —  Va 


(2-7) 


(2-8) 


are  obtained. 

Choosing  state  variables  as: 

Xi  =  w 
X2  =  ia 

Xl   -  -  j    X1      +  —  X2 

X2  =  -  5*  Xl   +  H  X2  +  1  Vfl 

L        L       L   a 
A  signal  flow  graph  of  the  system  is  given  in  Figure  2.5, 
with  1/s  indicating  an  integration  and  Xio  and  X20   as  initial 
conditions . 
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FIGURE  2.5:   Signal  flow  graph  for  linear  system 
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Using  Mason's  gain  rule  the  transfer  function  in  Laplace 
notation  is: 


G2(s) 


Kt 

JL 

Va 

2         RJ+Lf 
s      +   

JL 

s+Rf+KtKm 

JL 

(2-9) 


or  in  Bode  form 


Kt 


G  (s)  =  Rf+KrKm 


JL     2    RJ+Lf       . 

S    +  br-s-TTS  „      S   +  1 


Rf+K,K        Rf+K,K  -(0    1nx 

t  m  t  m  (2-10) 


Numerical  values  for  the  motor  parameters  and  load  were 
selected  from  Reference  1,  Chapter  7,  Problem  9. 
They  are : 

Armature  resistance  R  =  0.08  ohm 

Armature  inductance  L  =  0.008  Henry 

Motor  torgue  constant  K  =  0.95  Newton-meter/Amp 

Motor  e.m.f.  constant  K  =  0.95  Volt-sec/Rad 

m 


Equivalent  inertia  of  motor-load  J  =  1.38  Kg-m 


2 


Equivalent  viscous  friction  of   f  =  0.086  newton-meter 
motor-load  rad/sec 

Feedback  gain  K  =  0.318  Volt-sec/rad. 

By  using  these  values  and  computer  program  "LISA  360  A  PRO- 
GRAM FOR  LINEAR  SYSTEM  ANALYSIS  360  D-16  .  4  .  009",  the 
values  of  magnitude  in  "db"  and  phase  of  the  transfer  func- 
tion were  calculated.   From  these  data  Bode  plot  and  Nichols 
plot  were  obtained  for  limit  cycle  investigation.   Plots  were 
shown  in  Figures  2.6  and  2.7- 
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FIGURE  2.7:   Nichols  chart  for  linear  part 
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For  a  desired  value  of  speed  at  a  specified  load,  normal 
practice  is  to  compute  the  value  of  firing  angle  that  will 
produce  the  desired  output,  then  design  the  firing  circuit  to 
produce  a  firing  angle,  ot0  9    that  biases  the  rectifier  to  the 
desired  output.   Thus,  the  basic  design  provides  the  desired 
result  in  open  loop  operation.   The  feedback  loop  is  then 
designed  to  adjust  the  system  when  operating  conditions  change 
If  the  load 'is  changed  this  causes  a  change  in  speed  which 
causes  an  error,  and  the  error  is  amplified  and  changes  the 
firing  angle  by  an  amount  3  in  order  to  counteract  the  'change 
in  speed.   In  such  a  design  it  is  necessary  that  the  firing 
circuit  gain,  K  ,  be  chosen  to  give  the  desired  steady  state 
speed  regulation,  then  the  loop  must  be  checked  for  stability 
and  for  transient  response,  since  the  gain  is  set  on  the 
basis  of  steady  state  accuracy.   It  is  possible  to  use  an 
open  loop  approach  which  includes  the  nonlinear  characteris- 
tic of  the  rectifier. 

Prom  Figure  2.2, 

fi(s)  =  Gi (s)  N  G2(s) 
and 

B(s)  =  G2(s)  V&(s)  (2-11) 

Assuming  armature  voltage  is  a  constant,  that  is  the  average 
value  of  the  rectifier  output,  V   .   Using  equation  2-3  and 
motor  parameters  equation  2-11  becomes 

86  Vav 


fi(s)    s(sz  +  10.06s  +  82.37) 


2k 


By  applying  the  final  value  theorem  which  is 

lim  w(t)  =  lim  s^(s) 
t-*-°°        s+o 

The  steady  state  value  of  speed  was  obtained  as  below: 

u>   =  1.045V 

ss         av 

It  is  known  that  the  average  value  of  the  rectifier  output  is 
a  function  of  the  error  signal,  a,  and  the  sensitivity  of  the 
firing  circuit  to  the  error.   If  the  sensitivity  of  the  fir- 
ing circuit  is  improved  by  using  an  amplifier  and  compensation 
network,  better  results  can  be  obtained. 

As  mentioned  before,  operation  of  the  system  depends  on 
existence  of  a  difference  between  reference  speed  and  measured 
speed.   If  the  "control  gain"  of  the  system  is  defined  as  the 
tachometer  output  per  volt  applied  to  the  amplifier  when  the 
feedback  path  is  broken. 

Control  gain 

K,   u 
G  =  .  t   m 


a 


Now  when  the  loop  is  closed 


Then, 


a  =  V   -  K.  03 
r    t  m 


K,  a) 
G  =    t      m 


Vr  "  Vm 


K,  (o   =  V    G 
t  m    r 


1  +  G 


Thus,  the  greater  the  control  gain,  the  closer  will  be  K,uj 
to  V  .   For  high  accuracy,  high  gain  is  required. 
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III.   PREDICTION  OF  LIMIT  CYCLE  CHARACTERISTIC 

Procedures  for  analysis  of  the  linear  portion  of  the 
system  have  been  previously  established  by  using  the  trans- 
fer function  concept.   It  is  desirable  to  extend  these 
techniques  to  the  analysis  of  the  nonlinear  part  of  the  sys- 
tem.  Actually,  for  the  nonlinear  portion  there  is  generally 
no  exact  method  for  relating  the  output  to  the  input,  in 
closed  form.   Most  techniques  of  nonlinear  analysis  use  some 
sort  of  approximation,  that  is,  assumed  linearity  of  the 
system  over  the  desired  range  of  operation  to  be  considered. 

The  extension  of  the  transfer  function  representation 
is  more  difficult.   It  requires  that  the  nonlinearity  be 
described  by  an  equation  in  the  frequency  domain,  and  it  also 
requires  that,  this  equation  be  compatible  for  use  with  the 
transfer  functions  of  the  linear  components.   Practical  solu- 
tion to  these  requirements  is  the  describing  function  technique 
which  makes  the  approximation  in  the  frequency  domain. 

The  basic  philosophy  of  the  describing  function  is  that 
the  open-loop  control  system  is  excited  at  its  input  by  a 
sinusoidal  forcing  function.   The  output  of  the  nonlinearity 
is  periodic  and  can  be  represented  by  its  Fourier  series. 
This  is  illustrated  in  Figure  3.1. 

Certain  assumptions  must  be  made  in  the  process  of  this 
linearization,  that  the  input  to  the  nonlinear  component  is 
a  pure  sine  wave,  and  linear  portion  of  the  system  following 
the  nonlinearity  is  a  low-pass  filter,  which  is  commonly  true 
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in  most  physical  feedback  control  systems.   Thus  it  is  assumed 
that  the  only  terms  of  the  Fourier  series  that  are  passed  to 
the  output  of  the  system  are  the  D-C  term  and  the  fundamental 
frequency  term.   The  fundamental  frequency  term  in  the  Fourier 
series  has  the  same  frequency  as  the  input  signal  but  may  dif- 
fer in  amplitude  and  phase.   When  the  loop  is  closed  this 
harmonic  will  propagate  itself  around  the  loop  and  result  in 
a  continuous  oscillation  of  the  system.   Then  the  describing 
function  may  be  defined  as  the  ratio  of  the  fundamental  term 
in  the  Fourier  series  for  the  output  wave  to  the  input  sinu- 
soid considering  both  magnitude  and  phase. 

In  order  to  get  the  describing  function  for  a  nonlinear- 
ity,  the  output  of  this  nonlinearity  is  represented  by  its 
Fourier  series. 

A 
f(t)  =  ■—-  +   AiCoswt  +  BxSinwt  + +  A  Cosntot  +  B  Sinncot 

(3-D 
The  coefficients  of  the  series  are  given  by 


An   =  ~f~y 

4) 


T 
B  =  -Tjr- I         f(t)    Sinncot   dt  (3-2b) 

Jo 

The  definition  of  the  describing  function  requires  that  A 
and  B   for  n  >  2  are  negligible.   Thus  the  coefficients  of 
interest  are  Ai  and  Bi ,  and  they  define  the  fundamental  fre- 
quency output  as 

0(t)  =  |  A!2  +  Ba  2  |  ^   /  tarT^Ai/Bi)  (3-3) 
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Then  the  describing  function  for  the  system  as  illustrated 
in  Figure  3.1  is 

I  A!  2   +  B^l^ 

N  = /tan"1(Ai/Bi)  -  <j> 

B  (3-6) 

The  above  criteria  will  now  be  applied  to  the  system 
under  consideration. 

A  mathematical  description  of  the  output  wave  shape  of 
the  nonlinear  block  must  first  be  obtained  for  derivation  of 
the  describing  function.   After  that  the  Fourier  analysis  is 
worked  out  to  get  sinusoidal  components  of  it. 

Referring  to  Figure  3-2,  the  output  of  the  system  has  a 
certain  constant  speed,  and  a  ripple  which  is  superimposed  on 
this  constant  level.   The  first  harmonic  of  the  ripple  is  de- 
fined as  BSinwt,  and  the  difference  between  the  constant  level 
and  the  reference  speed  as 

ao  =  VR  "  VDC  " 

The  system  is  presumed  to  operate  as  follows : 

a)  When  the  constant  level  of  the  speed  plus  the  ripple 
is  less  than  the  reference  level  the  rectifier  bridge  is  "ON" 
between  time  a  to  b,  and  the  output  of  the  nonlinear  portion 
of  the  system  is  a  three-phase,  full-wave  rectified  voltage 
of  mean  value  A. 

b)  At  time  b  where  the  constant  speed  plus  the  ripple 
equals  the  reference  speed,  the  rectifier  bridge  receives  an 
"OFF"  signal. 
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BSincot 


f  (I:)  =ZAnCosnojt  +  ZBnSinnojt 


FIGURE    3.1. 


f(t) 


*»-  time 


REF 


FIGURE    3.2. 
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c)  The  rectifier  is  turned  off  at  time  equal  to  c,  a 
delay  of  1/240  second  after  time  b. 

d)  The  rectifier  remains  off  until  the  output  speed 
again  falls  below  the  reference.   (Time  c  to  T.)   Then  the 
output  of  the  nonlinear  block  may  be  summarized  as  follows 

f(t)  =  0         0<t<a 
=  A         a<t<b 

=  A   c"^    b<t<c 
c-b 

=  0  c<t<T 

The  times  a,  b5  and  c  can  be  defined  in  terms  of  the 
first  harmonic  of  the  ripple. 


b    = 

Sin" 

-1      a 

B 

03 

7T 

a  = 

-   b 

0) 

r»     = 

h     +    - 

1 

240 

To  find  the  coefficients  of  the  Fourier  series  equations 
3-2a,b  are  used.   The  integration  must  be  evaluated  over  one 
full  cycle  of  the  disturbance  frequency,  T. 

Procedure  is  as  follows : 


A   =  -|-  /   f  (t)  Cosnwt  dt 
n    T  / 
/  o 


^£  =  240b  +  1  -  240  t 
c-b 
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Then, 


A      =   M 

1  T 


J       Costot    dt+240b  j       Coscot    dt-240  /       t    Cosoot    dt 


Integrating   and   substituting   the    limits 


A      =   ^ 

1         Too 


SinioC-Sinaja+240bSin(oc-240bSinojb    -^-^  Cosooc 

0) 


-240cSintoc  +  -^-  Cosoob  +  240bSin(ob 
to 


Cancelling   terms,    the    above    reduces    to: 


A      =   M 

Al         Toj 


-    Sintob  +   (Cosoob-Cosooc ) 


Using  trigonomatric  identities  and  substituting  the  follow- 
ing values, 


T  a 

=  7TT-  ,  Sintob  =  q  and  Cosoob  = 


2tt 


B  . 


/B7^T7 
B 


The    final    form   of   the   A.    term   is    found   as    equation    3-5- 


A      =   -^ 
Al        ttB 


■  a   + 


240 


oo 


/ 


2      „,  2 


Bz-a 


1-Cos 


CO 


2^0 


+aSin 


oo 


2"4~0 


(3-5) 


Proceeding   in   a   like   manner   B.,    is    found   as    equation    3-6. 


Bl    "   ^B 


vC7 


-     B 


2         240 
a      +   


w 


a 


1-Cos 


^ 

240 

4 

-v . 

2    oin    0iir, 
-    a                240 

(3-6) 
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The  describing;  function  is  then 


^Ai2  +  Bx2 

N  =  


/e 


B  (3-7) 

where 

G  =  Tan"1  (A i /BO 

Equations  3-5  >  3-6,  and  3-7  are  suitable  for  evaluating 
the  numerical  value  of  the  describing  function.   Due  to  their 
complex  nature,  the  digital  computer  was  used  to  overcome  the 
labor  required  to  calculate  more  than  a  limited  number  of 
points . 

The  computer  program  with  numerical  data  is  contained  in 
Appendix  A. 

Since  the  describing  function,  which  is  derived  above  is 
dependent  upon  the  frequency,  the  magnitude  of  the  input,  B, 
and  the  difference  between  constant  level  of  the  output  and 
the  reference  speed  a0 ,  separate  sets  of  calculations  are 
necessary  for  each  choice  of  an  independent  variable.   For 
convenience  the  magnitude  of  the  input  was  chosen  as  the  in- 
dependent variable.   The  flow  chart  for  the  computer  program 
is  shown  in  Figure  3 -3> 

To  make  the  prediction  of  limit  cycle  characteristic  the 
describing  functions  were  drawn  on  the  Nichols  chart.   They 
are  shown  in  Figures  3-^  through  3*6. 

The  values  of  alpha  and  frequencies  for  evaluation  of 
describing  functions  were  obtained  from  the  simulation  of  the 
system  which  will  be  discussed  later. 
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FIND 
DELTQ     S 
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1 
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N  ,  -  -if  (db)  ,  3  ,    Angle 

-1 

' 

-» 

B  =  8  ■*  DELTAS 

Print    Output 
&  >  Nt  -  i/n  ,     Angle 


NO 


FIGURE  3.3 


Flow  chart  to  compute  value  of  describing 
function. 
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# 

Alpha 
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1 

2 
3 

0.01 

0.0055 

0.0035 

51.00  HZ 
56.61  HZ 
62.00  HZ 

0.0109 


_V5' 


-50 


db. 


0.00^2* 

(Phase) 


-l?D  -MS 


FIGURE    3.4: 


Describing  function  (-1/N)  for  limit  cycling 
thyristor  rectifier  bridge. 
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# 

Alpha 

Frequency 

1 

2 
3 

0.01 

0.0065 

0.0045 

45.00  HZ 
50.15  HZ 
55.00  HZ 

+165 


0.0D15        \6-0l2l 


FIGURE  3.5 


Describing  function  (-1/N)  for  limit  cycling 
thyristor  rectifier  bridge. 
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I 

115' 


-to 
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/go' 


-/75°  ~/?0' 

(phgse) 


/6S' 


FIGURE  3.6: 


Describing  function  (-1/N)  for  limit  cycling 
thyristor  rectifier  bridge. 
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These  values  were  tabulated  in  Table  3-1. 


REFERENCE    MEAN  VALUE 
r ^m IN  RAD/SEC   OF  OUTPUT    ALPHA    FREQUENCY 

31v     930  rpm     97 • 4882     97 • 4843    0.0055    56.61HZ 


32. 8V   984  rpm    103-1443    103.1378    O.OO65    50.15HZ 
33. 5V   1005  rpm    105.3459    105.337     0.0089    39.177HZ 

TABLE  3-1.   Tabulation  of  simulation  results  for  V   equal 
31,32.8  and  33-5  Volts .  V 

Having  derived  the  transfer  function  for  the  linear  por- 
tion of  the  system  and  describing  function  for  the  nonlinear 
part  of  the  system  there  are  now  a  number  of  ways  to  deter- 
mine the  condition  necessary  for  limit  cycling. 

The  characteristic  equation  for  the  system  is  written 
from  equation  2-3  as: 

1  +  Gi (s)  G2(s)  H(s)  N  =  0 
This  can  be  rearranged  as : 

G1(s)G2(s)H(s)  =  -  I  =  G(s) 

Thus  the  intersection  of  a  graphical  representation  of  G  with 
-1/N  determines  the  critical  point  at  which  limit  cycling  will 
occur.   Either  a  Nyquist  or  a  Nichols  plot  is  convenient  for 
this  purpose.   The  choice  of  the  Nichols  plot  for  these  ana- 
lyses is  arbitrary. 
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To  determine  the  location  of  the  limit  cycle  the  des- 
cribing function  curves  were  superimposed  on  the  frequency 
response  curve  of  the  linear  part  for  frequency  39.177  HZ 
with  alpha  0.0089,  for  frequency  50.15  HZ  with  alpha  0.0065, 
and  for  frequency  56. 6l  HZ  with  alpha  0.0055.   The  graphs 
are  illustrated  in  Figures  3-7  through  3.9. 

The  values  of  frequency  and  amplitude  were  measured  at 
intersections  of  G  with  -1/N.  The  results  are  tabulated  in 
Table  3-2. 


V 


Amplitude  of   Frequency  of     Frequency  of 
r -1/N -1/N G 

31.  0.0063  56.61  HZ  56.65  HZ 


32.8  0.0072  50.15  HZ  50.13  HZ 


33.5  0.0095  39.17  HZ  39.23  HZ 


TABLE  3-2.   Tabulation  of  the  observed  limit  cycle  frequency 
and  amplitude. 


From  the  comparison  of  the  results  in  Table  3-2,  it  can 
be  seen  that  the  frequency  of  the  -1/N  curve  and  frequency  of 
the  G  curve  are  very  close  to  each  other.   This  is  true  for 
all  three  cases. 

These  results  were  obtained  when  alpha  was  equal  to  the 
values  that  are  indicated  in  Table  3-1-   Corresponding  to 
these  alpha  values  the  amplitude  of  the  limit  cycle  was  found 
as  indicated  in  Table  3-2. 
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f=  II.3H3 
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FIGURE    3.7. 
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f  =  n.3  Ht. 


VR  =  32.8  Volts 


Alpha 

0.01 

0.0065 

0.0045 


Frequency 

45.00  HZ 
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165° 


FIGURE    3.8 
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VR  =  33.5  Volts 

Alpha    Frequency 


1 
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35.00  HZ 
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US 


FIGURE  3.9. 
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Thus  there  is  good  agreement  with  simulation  results  in 
both  frequencies  of  the  limit  cycle  and  the  values  of  alpha 
predicted  by  the  describing  function.   This  indicates  the 
validity  of  the  describing  function  in  prediction  -of  the 
limit  cycle  of  the  system. 
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IV.   SIMULATION  OF  THE  SYSTEM  AND  SUBHARMONIC  EFFECTS 

To  determine  the  validity  of  the  describing  function 
and  to  see  the  performance  of  the  system  in  the  time  domain, 
the  system  was  simulated  on  the  digital  computer.   The  flow 
chart  for  the  simulation  is  shown  in  Figure  4.1.   Naval 
Postgraduate  School  computer  library  subroutine  "INTEG  2" 
was  used  to  solve  the  differential  equations  required  in  the 
simulation.   The  complete  program  is  contained  in  Appendix  B. 

The  supply  voltage  was  considered  to  be  a  three-phase 
60  cycle  A-C  source  with  a  peak  amplitude  of  115  volts  through- 
out the  simulation. 

The  forcing  voltage  for  the  system,  that  is  applied  to 
the  armature  of  the  D-C  motor  is  a  three-phase,  full-wave, 
rectified  voltage  during  time  a  to  b,  a  portion  of  a  sine 
wave  for  the  appropriate  conduction  phase  during  time  b  to 
c,  and  zero  voltage  with  no  reverse  current  through  the  in- 
ductance during  c  to  T  of  Figure  3.2. 

Graphical  output  from  this  program  are  shown  in  Figure  4.2 
illustrating  the  input  voltage  wave  form,  and  Figures  4.3 
through  4.7  illustrating  the  output  speed  wave  forms. 

Simulation  was  performed  by  changing  the  reference  from 
930  r.p.m.  to  1260  r.p.m.   These  are  shown  in  Table  4-1. 

When  a  linear  system  is  driven  with  a  sinusoidal  input 
it  is  known  that  the  steady  state  output  is  also  a  sinusoid, 
differing  from  the  input  in  amplitude  and  phase  but  having 
precisely  the  same  frequency.   For  linear  systems,  that  is 
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FIGURE  4.1.   Flow  chart  for  simulation 
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FIGURE  4-2.   Simulation  results,  input  voltage. 
Reference  speed  1050  rpm. 
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FIGURE  4-3.   Simulation  results,  output  speed. 
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"FIGURE  4-4.   Simulation  results,  output  speed 
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FIGURE  4-5.   Simulation  results,  output  speed. 
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FIGURE  4-6.   Simulation  results,  output  speed. 
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FIGURE  4-7.   Simulation  results t    output  speed. 
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REFERENCE 
VOLTAGE 

SPEED 

PERIOD  OF 
LIMIT  CYCLE 

FREQUENCY  OF 
LIMIT  CYCLE 

DOES  SUB- 
HARMONIC 
EXIST? 

31 

volt 

930 

rpm 

0.017663 

Sec 

56.61 

HZ 

YES 

32 

volt 

960 

rpm 

0.017631 

Sec 

56.71 

HZ 

YES 

32.3 

volt 

969 

rpm 

0.021173 

Sec 

47.28 

HZ 

YES 

32.5 

volt 

975 

rpm 

0.019924 

Sec 

50.14 

HZ 

YES 

32.6 

volt 

978 

rpm 

0.0195 

Sec 

51.28 

HZ 

YES 

32.8 

volt 

984 

rpm 

0.01994 

Sec 

50.13 

HZ 

YES 

33. 

volt 

990 

rpm 

0.02391 

Sec 

41.82 

HZ 

YES 

33-3 

volt 

999 

rpm 

0.0223 

Sec 

44.576 

HZ 

YES 

33.5 

volt 

1005 

rpm 

0.025525 

Sec 

39.177 

HZ 

YES 

33.8 

volt 

1014 

rpm 

0.02545 

Sec 

39.293 

HZ 

YES 

3^. 

volt 

1020 

rpm 

0.027838 

Sec 

35.923 

HZ 

YES 

35. 

volt 

1050 

rpm 

0.0394 

Sec 

25.4 

HZ 

YES 

37. 

volt 

1110 

rpm- 

NO 

LIMIT 

CYCLE 

40. 

volt 

1200 

rpm 

NO 

LIMIT 

CYCLE 

42. 

volt 

1260 

rpm 

NO 

LIMIT 

CYCLE 

TABLE  4-1.   Tabulation  of  simulation  results. 
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true  at  all  frequencies.   In  the  case  of  a  nonlinear  system 
a  sinusoidal  input  does  not  guarantee  an  output  which  is  a 
sinusoid.   In  the  usual  case,  the  output  wave  contains  har- 
monic frequency  components  which  are  integral  multiples  of 
the  forcing  frequency.   In  some  cases,  the  output  wave  may 
contain  one  or  more  frequencies  which  are  lower  than  the  for- 
cing frequency.   Characteristically,  these  frequencies  are 
integral  submultiples  of  the  forcing  frequency  which  are 
called  subharmonic  frequencies.   The  amplitude  of  the  sub- 
harmonic  component  may  be  small  compared  with  the  amplitude 
of  the  component  at  the  driving  frequency,  or  it  may  be  so 
large  that  the  driving  frequency  may  be  neglected  for  partic- 
ular computations. 

If  the  output  wave  shapes  in  Figure  4.3  through  4.7  are 
inspected,  it  is  observed  that  the  wave  forms  have  a  ripple 
instability  which  is  small  in  amplitude  compared  to  the  ave- 
rage ripple  amplitude.   One  can  think  this  instability  is 
caused  by  subharmonic  effects  of  the  system.   It  can  be  ob- 
served by  inspection,  but  it  is  not  possible  to  obtain  its 
frequency  and  its  amplitude  from  the  output  wave  shapes,  be- 
cause they  do  not  have  a  regular  pattern.   To  figure  out  this 
instability,  the  Fourier  analysis  of  the  output  wave  shapes 
for  reference  speed  930,  984,  and  1005  r.p.m.  were  made. 
During  the  Fourier  analysis,  several  limit  cycle  periods  were 
taken  as  a  fundamental  period  of  the  Fourier,  then  the  same 
harmonics  of  these  fundamental  frequencies  which  are  also  sub- 
harmonics  of  the  limit  cycle  were  obtained.   From  these  re- 
sults, amplitude  and  the  order  of  subharmonic  components  of 
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the  limit  cycle  in  the  output  wave  were  calculated.   These 
were  tabulated  in  Table  4-2  through  4-4. 

In  Figure  4.8  through  4.10,  the  amplitude  of  the  sub- 
harmonics  was  drawn  to  see  the  relative  heights  of  the  com- 
ponents.  In  Figure  4.8,  it  is  seen  that,  at  reference  speed 
930  r.p.m.  1/8  and  1/16  subharmonics  have  larger  amplitude 
compared  to  others.   Their  amplitudes  are  approximately  15.6 
percent  of  the  output  ripple  amplitude.   In  Figure  4.9  which 
was  drawn  for  reference  speed  984  r.p.m.,  1/2  subharmonic 
has  larger  amplitude  and  it  is  equal  25-3  percent  of  the  out- 
put ripple.   In  Figure  4.10,  which  is  for  reference  speed 
1005  r.p.m.  the  largest  subharmonic  is  1/7.   It  has  amplitude 
16.1  percent  of  the  ripple  amplitude. 
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V.   CONCLUSION  AND  RECOMMENDATION  FOR  FURTHER  STUDY 

The  main  conclusion  that  might  be  reached  in  reviewing 
the  experimental  results  is  that  the  D-C  motor  at  constant 
field  current  operates  at  an  almost  constant  speed,  and  its 
speed  may  be  controlled  over  wide  ranges  by  armature-voltage 
control.   This  can  be  seen  from  Figures  4.3  through  4.7.   If 
these  figures  are  inspected  closely  it  can  be  observed  that 
the  ripple  amplitudes  are  very  small  compared  to  the  output 
value.   They  are  approximately  within  the  ±  0.025  percent  of 
the  desired  value  of  the  output.   This  may  be  acceptable  for 
most  of  the  applications  in  the  industry. 

Armature-voltage  control  was  performed  by  adjusting  the 
conduction  phase  of  the  thyristor-rectif ier  bridge.   This 
causes  the  limit  cycling  in  the  control  loop. 

The  describing  function  derived  in  Chapter  III  provides 
a  mean  by  which  the  nonlinearity  of  the  system  can  be  repre- 
sented and  the  overall  performance  of  the  system  analyzed. 
It  is  possible  to  predict  the  approximate  frequency  and  ampli- 
tude of  the  limit  cycle  as  illustrated  in  Figures  3«7  through 
3.9«   In  order  to  predict  the  response  more  precisely,  it  is 
necessary  to  predetermine  the  value  of  alpha,  that  might  be 
done  by  simulating  the  system  and  measuring  the  output  mean 
value  then  taking  the  difference  between  output  mean  value  and 
reference. 

One  tentative  conclusion  that  can  be  achieved  by  examina- 
tion of  the  experimental  data  is  that,  under  certain  conditions 
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the  f orced-limit-cycling  speed  control  system  will  have  sub- 
harmonic  ripple  instability. 

The  describing  function  derived  accounts  only  for  the 
first  harmonic  of  the  output  voltage  from  the  nonlinear  por- 
tion of  the  system  and  assumes  that  higher  harmonics  are 
negligible.   However,  as  seen  in  Figures  4.3  through  4.7  the 
other  harmonics  can  exist.   In  some  cases  these  harmonics 
can  be  submultiples  of  the  limit  cycle  frequency. 

In  this  study  Fourier  analysis  was  performed  to  see 
whether  the  subharmonic  exists  or  not.   As  a  result  Tables 
4-2  through  4-4  and  Figures  4.8  through  4.10  were  obtained. 
These  show  that  the  ripple  instability  is  caused  by  subhar- 
monic effects  of  the  system  nonlinearity .   But  this  is  not 
good  enough  to  predict  the  existence  of  the  subharmonic  in 
the  output  of  the  system.   It  can  give  some  insight  about 
the  nature  of  the  problem.   There  are  no  general  rules  avail- 
able which  define  the  conditions  necessary  for  the  occurrence 
of  the  subharmonic  oscillation.   For  mathematical  analysis 
several  of  the  references  noted  in  the  List  of  References 
suggest  that  a  possible  method  of  prediction  of  the  subhar- 
monic ripple  instability  in  the  system  would  be  to  employ  a 
Dual  Input  Describing  Function.   Any  further  study  should 
give  some  consideration  to  these  analytical  methods  to  under- 
stand ripple  instability  and  existence  of  the  subharmonic. 

Another  area  for  study  is  to  increase  the  gain  and  see 
the  effect  of  the  gain  on  the  ripple  instability.   This  may 
require  some  sort  of  compensation  network  to  eliminate  the 
instability  in  the  output  for  a  range  of  gain  values. 
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APPENDIX  A 

Computer  program  for  computing  value  of  describing 
curves.   Numerical  data  for  -1/N  curves  plotted  in  Figure  3-^ 
through  3«9- 
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DESCRIBING  FUNCTION  DATA 
FREQUENCY  =  51.000 

ALPHA  =  0.0100 


B 

N 

PHASE 

-1/N(DB) 

0.0101 

1401.8879 

29.27 

-62.9343 

0.0102 

106G.7751 

32.61 

-60.5125 

0.0102 

805.8579 

36.95 

-58.  1252 

0.01C3 

600.5212 

43.26 

-55.5706 

0.0104 

433. 9397 

53.37 

-5  2.7^86 

0.0105 

308,9268 

70.8  7 

-49.7971 

0.0106 

244-2780 

99,49 

-47,7577 

0,0106 

257.4648 

130.96 

-48.2144 

0,0107 

324-C369 

152. 10 

-50-2119 

0*0108 

410-26*9 

164-16 

-52,2613 

o.oiov 

50U.6140 

171.37 

-53.9901 

0.0110 

589.5430 

176.04 

-55.4103 

0.0110 

675. 1055 

179.29 

-56.5874 

0.0111 

756.6646 

-178.34 

-5  7.5781 

0.0112 

834.0947 

-176.53 

-58.4243 

0.0113 

907.4871 

-175.10 

-59.1568 

0.0114 

977.0127 

-17  3.95 

-5  9.79  80 

0.0114 

1042.8  850 

-17  2.99 

-60.3647 

0.0115 

1105.3154 

-172. 19 

-60.8697 

0.0116 

116'+. 5139 

-171.51 

-61.3229 

0.0117 

1220.6768 

-170.92 

-61.7320 

0.0118 

1273.9927 

-170.41 

-62. 1033 

0.0118 

1324.6333 

-169.95 

-62.4419 

0.0119 

1372.7588 

-169.55 

-62.7519 

0.0120 

1418.5159 

-169.19 

-63.0367 

0.0121 

1402.0398 

-168.67 

-63.2992 

0.0122 

1503.4609 

-168.57 

-63.5418 

0.0122 

1542.8931 

-168.31 

-63.7667 

0.0123 

1580.4492 

-168.06 

-6  3.97  56 

0.0124 

1616.2246 

-167.84 

-6^.1700 

0.0125 

1650.3162 

-167.63 

-64.3513 

0.0126 

1682.8105 

-167.44 

-64.5207 

0.0126 

1713.7903 

-167.26 

-64.6791 

0.0127 

1743.3315 

-167.10 

-64.82 16 

0.0128 

1771.5032 

-166.94 

-64.9668 

0.0129 

1798.3711 

-166.80 

-65.0976 

0.0130 

1824.0007 

-166.66 

-65.2205 

0.0130 

1848.4512 

-166.53   • 

-65.3362 

0.0131 

1871.7751 

-166.41 

-6  5.4451 

0.0132 

1894.0234 

-166.30 

-65.5477 

0.0133 

1915.2488 

-166.  19 

-6  5.6445 

0.0134 

1935.4956 

-166.09 

-65.7358 

0.0134 

1954.8059 

-165.99 

-65.8221 

0.0135 

1973.2205 

-165.90 

-65.9035 

0.0136 

1990.7820 

-165.81 

-65.9805 

0.0137 

2007.5239 

-165.73 

-66.05  32 

0.0138 

2023.4832 

-165.65 

-66.1220 

0.0138 

2038.6892 

-165.57 

-66.1870 

0.0139 

2053. 1780 

-165.50 

-66.2485 

0.0140 

2066.9778 

-165.43 

-66.3067 

DESCRIBING  FUNCTION  DATA 
FREQUENCY  =  56.610 


ALPHA  =  0.0055 


N 


0.0056 

0.0056 

0.0057 

0.0057 

0.0058 

0.0058 

0.0059 

0.0059 

0.0060 

0.0060 

0.0061 

0.0061 

0.0062 

0.0063 

0.0063 

0.0064 

0.0064 

0.0065 

0.0065 

0.0066 

0.0066 

0.0067 

0.0067 

0.0068 

0.0069 

0.0069 

0.0070 

0.0070 

0.0071 

0.0071 

0.0072 

0.0072 

0.0073 

0.0073 

0.0074 

0.0074 

0.0075 

0.0076 

0.0076 

0.0077 

0.0077 

0.0078 

0.0078 

0.0079 

0.0079 

0.0080 

0.0080 

0.0081 

0.0081 

0.0082 


2811 
2136 
1634. 
1233. 
911. 
671. 
542. 
547. 
652. 
802. 
965. 
1127. 
12  84. 
1434. 
1576. 
1710. 
1837. 
1956, 
2069, 
2175 
2276 
2370 
2460 
2544 
2624 
2700 
2771 
2839 
2903 
2963 
3021 
3075 
3126 
3175 
3221 
3264 
3306 
3345 
3382 
3417 
3450 
3481 
3510 
3538 
3564 
3589 
3613 
3635 
3656 
3675 


.1938 
.1204 
.  1643 
.0056 
.0759 
.7336 
2371 
1372 
74  76 
6943 
2810 
5667 
5220 
2671 
2329 
,4219 
,0771 
.5686 
.2861 
,6387 
.0090 
.7688 
.2668 
.8293 
.7527 
.3184 
.7900 
.3950 
.3601 
.9014 
.2019 
.4402 
.7830 
.3940 
.4092 
.9675 
.  1929 
.2063 
.1211 
.0393 
.0581 
.2764 
.7781 
.6392 
.9492 
.7737 
.1768 
.2361 
.0046 
.5469 


PHASE 


32.44 
36.08 
40.74 
'+7.36 
57.52 
74.02 
99.22 
127.39 
148.35 
161.28 
169.32 
174.63 
17  8.34 
-178.93 
-176.85 
-175.20 
-173.87 
-172.78 
-171.85 
-171.07 
-170.39 
-169.80 
-169.27 
-168.81 
-168.40 
-168.02 
-167.69 
-167.38 
-167. 
-166, 
-166 
-166.38 
-166.18 
-165.99 
-165.81 
-165.64 
-165.49 
-165.34 
-165.20 
-165.07 
-164.95 
-164.83 
-164.72 
-164.61 
-164.51 
-164.42 
-164.33 
-164.24 
-164.16 
-164.08 


10 
84 
60 


-1/N(DB) 


-68.9778 
-66.5925 
-64.2659 

-61.8193 

-59.1911 

-56.5^39 

-54.6838 

-54.7619 

-56.2949 

-58.0910 

-59.6931 

-61.0428 

-62.1748 

-63.1326 

-63.9524 

-64.6620 

-65.282  5 

-65.8299 

-66.3164 

-66.7517 

-67.1435 

-67.4978 

-67.8196 

-68.1132 

-68.3818 

-68.6283 

-68.8552 

-69.0645 

-69.2580 

-69.4373 

-69.6036 

-69.7581 

-69.9019 

-70.0359 

-70.1609 

-70.2776 

-70.3866 

-70.48  84 

-70.5838 

-70.6730 

-70.7565 

-70.8348 

-70.9081 

-70.9767 

-71.0411 

-71.1013 

-71.1578 

-71.2106 

-71.2601 

-71.3064 


65 


DESCRIBING  FUNCTION  DATA 
FREQUENCY  =  62.000 

ALPHA  =  0.0035 


B 

N 

PHASE 

-1/N(DB) 

0.0035 

^775.6094 

35.54 

-73.5806 

0.00  36 

3626.2739 

39.57 

-71.  1892 

0.0036 

2777.4131 

44.70 

-68.8728 

0.0037 

2106.3257 

51.87 

-66.4705 

0.0037 

1577.0222 

62.62 

-63.9568 

0.003  8 

1194.6785 

79.28 

-61.5450 

0.0038 

99o.3599 

103.00 

-59.9683 

0.0038 

1004.8066 

128.44 

-60.0417 

U.0039 

1165. 1162 

147.89 

-61.3274 

0.0039 

1396.2275 

160.52 

-62.8991 

0.0040 

1649. ol23 

168.68 

-64. 3476 

0.0040 

1903.7549 

174.20 

-65.5922 

0.0040 

2149.7078 

178.11 

-66.6476 

0.0041 

2383.9434 

-178.98 

-67.5459 

0.0041 

2605.2627 

-176.74 

-68.3170 

0.0042 

2813.5598 

-174.96 

-63.9351 

0.0042 

3009. 1775 

-173.52 

-69.5690 

0.004  3 

3192.7107 

-172.32 

-70.0332 

0.0043 

3364.8289 

-171.32 

-70.5393 

0.0043 

3526.2234 

-170.46 

-70.9462 

0.00^4 

3677.5554 

-169.72 

-71.3112 

0.0044 

3819.4849 

-169.07 

-71.6401 

0.0045 

3952.6213 

-168.50 

-71.9377 

0.0045 

4077.5251 

-167.99 

-72.2079 

0.0045 

4194.7266 

-167.53 

-72.4541 

0.0046 

430^.7109 

-167. 12 

-72.6789 

0.0046 

4407.9453 

-166.75 

-72.8347 

0.0047 

4504.8438 

-166.41 

-73.0736 

0.0047 

4595.7930 

-166.10 

-73.2472 

0.0048 

4681.  1602 

-165.82 

-73.4071 

0.0048 

4761.2773 

-165.56 

-73.5545 

0.0046 

4636.4609 

-165.31 

-73.6905 

0.0049 

4906.9883 

-165.09 

-73.8163 

0.0049 

4973. 1328 

-164.88 

-73.9326 

0.0050 

5035. 1602 

-164.68 

-74.0403 

0.0050 

5093.2969 

-164.50 

-74.1400 

0.0050 

5147.7500 

.  -164.33 

-74.232  3 

0.0051 

5196.7344 

-164.17 

-74.3179 

0.0051 

5246.4414 

-164.01 

-74.3973 

0.0052 

5291.0391 

-163.67 

-74.47  08 

0.0052 

5332.7109 

-163.73 

-7^.5390 

0.0053 

5371.5938 

-163.60 

-74.60  21 

0.0053 

5407.8516 

-163.48 

-74.6o05 

0.0053 

5441.6172 

-163.36 

-74.7146 

0.0054 

5473.0156 

-163.25 

-74.7o45 

0.0054 

5502. 1797 

-163. 15 

-74.8107 

0.0055 

5529.2188 

-163.05 

-74.8533 

0.0055 

5554.2422 

-162.95 

-74.89  2  5 

0.0056 

5577.3398 

-162.86 

-74.9285 

0.0056 

5598.6250 

-162.77 

-74.9616 

DESCRIBING  FUNCTION  DATA 
FREQUENCY  =  45.000 

ALPHA  =  0.0100 


B 

N 

PHASE 

-1/N(DB 

0.0101 

1233.4875 

25.95 

-61.8227 

0.0101 

922.4802 

29.06 

-59.2991 

0.0102 

689.1995 

33.25 

-56.7669 

0.0103 

500.3091 

39.63 

-53.9848 

0.0103 

346.5730 

50.63 

-50.7959 

0.0104 

233.9463 

71.83 

-47. 3823 

0.0105 

190.4*55 

108.20 

-45.5954 

0.0105 

230. 1004 

141.57 

-47.2383 

0.0106 

310.6777 

159.65 

-49.8462 

0.0106 

401.2224 

169.15 

-52.0677 

0.0107 

491.7175 

174.73 

-53.8343 

0.0108 

579.1682 

178.35 

-55.2561 

o.oioa 

662.6948 

-179.12 

-56.4263 

0.0109 

742. 1340 

-177.26 

-57.4096 

0.0110 

817.5820 

-175.84 

-53.2506 

0.0110 

889.2295 

-174.71 

-58.9803 

0.0111 

957.2971 

-173.79 

-59.6209 

0.0112 

1022.0098 

-173.03 

-60.  1891 

0.0112 

1083.5806 

-172.39 

-60.6972 

0.0113 

1142.2085 

-171.84 

-61.  1549 

0.0114 

119b. 0796 

-171.37 

-61.5697 

0.0114 

1251.3662 

-170.95 

-61.9477 

0.0115 

1302.2212 

-170.59 

-62.2937 

0.0116 

1350.7869 

-170.26 

-62.6117 

0.0116 

1397. 1980 

-169.97 

-62.9U52 

0.0117 

1441.5750 

-169.70 

-63.1767 

0.0117 

1484.0293 

-169.46 

-63.4288 

0.0118 

1524.6650 

-169.25 

-63.6635 

0.0119 

1563.5762 

-169.05 

-6  3.88  24 

0.0119 

1600.8506 

-168.86 

-64.0870 

0.0120 

1636.5732 

-168.69 

-64.2787 

0.0121 

1670.8147 

-168.54 

-64.45  86 

0.0121 

1703.6563 

-168.39 

-64.6276 

0.0122 

1735.1553 

-168.25 

-64.7863 

0.0123 

1765.3767 

-168.13 

-64.9368 

0.0123 

1794.3813 

-163.01 

-65.0763 

0.0124 

1822.2202 

-167.90 

-65.2120 

0.0125 

1848.9495 

-167.79 

-65.3385 

0.0125 

1874.6118 

-167.69 

-6  5.45  82 

0.0126 

1899.2578 

-167.60 

-65.5717 

0.0126 

1922.9282 

-167.51 

-65.6793 

0.0127 

1945.6631 

-167.43 

-65.7813 

0.0128 

1967.5032 

-167.35 

-65.8783 

0.0128 

1988.4817 

-167.27 

-65.9704 

0.0129 

2008.6379 

-167.20 

-66.0580 

0.0130 

2027.9993 

-167.13 

-66.1414 

0.0130 

2046.5977 

-167.06 

-66.2206 

0.0131 

2064.4663 

-167.00 

-66.2962 

0.0132 

2081.6311 

-166.94 

-66.36  81 

0.0132 

2098. 1150 

-166.88 

-66.4366 

DESCRIBING  FUNCTION  DATA 
FREQUENCY  =  50.150 

ALPHA  =  0.0065 


N 


PHASE 


-1/N(DB) 


0.0066 

0.0066 

0.0067 

0.0067 

0.0068 

0.0068 

0.0069 

0.0069 

0.0070 

0.0070 

0.0071 

0.0071 

0.0072 

0.0072 

0.0073 

0.0073 

0.0074 

0.0074 

0.0075 

0.0075 

0.0076 

0.0076 

0.0077 

0.0077 

0.0078 

0.0078 

0.0079 

0.0079 

0.0080 

0.0080 

0.0081 

0.0081 

0.0082 

0.0082 

0.0083 

0.0083 

0.0084 

0.0084 

0.0085 

0.0085 

0.0086 

0.0086 

0.0087 

0.0087 

0.0088 

0.0088 

0.0089 

0.0089 

0.0090 

0.0090 


2121.3572 
1603.3286 
1215.9641 
903.6431 
650.0144 
459.8135 
363.3149 
387.8416 
493.0381 
626.7251 
765.8857 
902.5076 
1033.8213 
1158.9558 
1277.7654 
1390.4126 
1497.1785 
1598.3855 
1694.3630 
1785.4329 
1871.8958 
1954.0325 
2032. 1028 
2106.3528 
2177.0103 
2244.2700 
2308.3335 
2369.3716 
2427.5530 
2483.0288 
2535.9387 
2586.4155 
2634.5789 
2680.5549 
2724.4355 
2766.3323 
2806.3376 
2844.5381 
2881.0166 
2915.8535 
2949.1211 
2980.8892 
3011.2258 
3040.1951 
3067.8445 
3094.2417 
3119.4341 
3143.4761 
3166.4092 
3188.2827 


28.79 
32.09 
36.40 
42.70 


86 
66 

,08 
02 


52, 
70 
100 
132, 
152.97 
164.77 
171.79 
176.33 
179.48 
-178.21 
-176.45 
-175.06 
-173.94 
-173.01 
-172.23 
-171.56 
-170.99 
-170.49 
-170.05 
-169.66 
-169.31 
-168.99 
-168.70 
-168.44 
-168.20 
-167.98 
-167.78 
-167.60 
-167.42 
-167.26 
-167.11 
-166.97 
-166.84 
-166.71 
-166.59 
-166.48 
-166.38 
-166.28 
-166.18 
-166.09 
-166.01 
-165.93 
-165.85 
-165.77 
-165.70 
-165.63 


-62 

-62 
-63 
-64 


66.5323 
64.1004 
61.6984 
59.1199 
56.2584 
53.2516 
51.2056 
51.7731 
53.6576 
5  5.9415 
•57.6833 
•59.1090- 
•60.2889 
■61.2813 
1290 
8629 
5054 
0736 
-64.5801 
-65.0349 
-65.4456 
-65.8186 
-66.1589 
-66.4706 
-66.7572 
-67.0215 
-67.2660 
-67.4926 
-67.7034 
-67.8996 
-68.0826 
-68.2540 
-68.4142 
-68.5645 
-68.7055 
-68.8381 
-68.9628 
-69.0802 
-69.1909 
-69.2953 
-69.3938 
-69.4869 
-69.5748 
-69.6580 
-69.7367 
-69.8111 
-69.8815 
-69.9482 
-70.0113 
-70.0711 


DESCRIBING  FUNCTION  DATA 
FREQUENCY  =  55.000 

ALPHA  =  0.0045 


B 

N 

PHASE 

-1/N(DB) 

0.0045 

3347.0078 

31.52 

-70.4931 

0.0046 

2541.6  218 

35.07 

-68.1022 

0.0046 

1941.8220 

39.62 

-65.7642 

0.0047 

1461.2092 

46.12 

-63.2942 

0.0047 

1074.0417 

56.20 

-60.6204 

0.C04& 

784.7676 

7  2.63 

-57.8948 

0.0048 

628.5144 

9  8.75 

-55.9663 

0.0048 

636.3621 

127.65 

-56.  1013 

0.0049 

771.3  501 

149.06 

-57.7450 

0.0049 

950.6079 

161.89 

-59.61^7 

0.0050 

1155.9346 

169.77 

-61.2587 

0.0050 

1354.2449 

174.94 

-62.6339 

0.0050 

1545.8279 

178.54 

-63.7332 

0.0051 

1728.6040 

-178.82 

-64.7539 

0.0051 

1901.9856 

-176.80 

-65.5841 

0.0052 

2  066.02  34 

-175.21 

-66.3027 

0.0052 

2221.0366 

-173.92 

-66.9311 

0.0053 

2367.4832 

-17  2.66 

-67.4857 

0.0053 

2505.6259 

-171.97 

-67.9790 

0.0053 

2636. 5632 

-171.21 

-68.4207 

0.0054 

2760.1492 

-170.56 

-68.8166 

0.00  54 

2677.0281 

-169.99 

-69.17  89 

0.0055 

2987.6121 

-169.43 

-69. 5065 

0.0055 

3092.2815 

-169.03 

-69.8056 

0.0055 

3191.3955 

-168.63 

-70.0796 

0.0056 

3265.2859 

-168.27 

-70.3315 

0.0056 

3374.2478 

-167.95 

-70.56  35 

0.00  57 

3456.5313 

-167.65 

-70.7779 

0.0057 

3536.5295 

-167.38 

-70.9764 

0.0058 

3614.3521 

-167.13 

-71.1606 

0.0058 

3686.2698 

-166.90 

-71.3317 

0.0058 

3754.4946 

-166.69 

-71.4910 

0.0059 

3819.2246 

-166.49 

-71.6395 

0.0059 

3860.6392 

-166.31 

-71.7761 

0.0060 

393d. 9158 

-16  6.14 

-71.9075 

0.0060 

3994.2109 

-165.98 

-72.0236 

0.0060 

4046.6799 

-165.83 

-72. 1420 

0.0061 

4096.4609 

-165.69 

-72.2482 

0.006  1 

4143.6336 

-165.55  • 

-72.3477 

0.0062 

4186.4727 

-165.43 

-  7  2 . 44 1 1 

0.0062 

4230.9531 

-165.31 

-72.5288 

0.0063 

4271.2305 

-165. 19 

-72.6111 

0.0063 

4309.3945 

-165.09 

-72.6883 

0.0063 

4345.5703 

-164.9  8 

-72.7609 

0.0064 

4379.8231 

-164.89 

-72.8291 

0.0064 

4412.2695 

-164.80 

-72.6932 

0.0065 

4442.9688 

-164.71 

-72.9535 

0.0065 

447  2.0  0  76 

-164.62 

-73.0100 

0.0065 

4499.4570 

-164.54 

-73.0632 

0.0066 

4525.3906 

—164.46 

-73.1131 

DESCRIBING  FUNCTION  DATA 
FREQUENCY  =  35.000 

ALPHA  =  0.0100 


B 

N 

PHASE 

-1/N(DB) 

0.0100 

923.5095 

20.57 

-59.3038 

0.0101 

657.7041 

23.60 

-56. 3606 

C.01C1 

-+57.7996 

23.24 

-53.2135 

0.0102 

296. 1536 

36.81 

-49.4303 

0.0102 

170.5451 

57.  14 

-44.6368 

0.0103 

118.4211 

107.9  0 

-41.46  86 

0.0103 

175. 1903 

151.02 

-44.8702 

0.0104 

267.0454 

167.  16 

-48.5317 

0.0104 

361.1838 

174.34 

-51.1546 

0.0105 

451.8193 

173.29 

-53.0993 

0.0105 

537.9329 

-179.23 

-54.6146 

0.0105 

619.5305 

-177.52 

-55.8412 

0.0106 

696.8745 

-176. 28 

-56.8631 

0.0106 

770.2817 

-175.33 

-57.7330 

0.0107 

840.0630 

-174.57 

-58.4862 

0.0107 

906.5017 

-173.97 

-59.  147  + 

0.0108 

969.8557 

-173.46 

-59.7341 

0.0103 

1030.3525 

-173.04 

-60.2597 

0.0109 

1086. 1973 

-172. 68 

-6  0.7341 

0.0109 

1143.5 71B 

-172.36 

-61  .1653 

0.0110 

1196. o338 

-172.09 

-61.5592 

0.0110 

1247.5388 

-171.84 

-61.921 1 

0.0110 

12  96.4126 

-171.63 

-62.25<v9 

0.0111 

1343.3779 

-171.43 

-62.5639 

0.0111 

1388.5435 

-171.26 

-62.8512 

0.0112 

1432.0044 

-171.10 

-63.  1169 

0.0112 

1473.6572 

-170.95 

-63. 3691 

0.0113 

1514.1843 

-170.82 

-6  3.60  36 

0.0113 

1553.0566 

-170.70 

-63.8237 

0.0114 

1590.5574 

-170.59 

-64.0310 

0.0114 

1626.7417 

-170.48 

-64.2264 

0.0114 

1661.6733 

-170.39 

-64.4109 

0.0115 

1695.4109 

-170.30 

-64. 5855 

0.0115 

1728.0044 

-170.21 

-64.7509 

0.0116 

1759.5068 

-170.13 

-64.90  7  8 

0.0116 

1769.9636 

-170.06 

-65.0569 

0.0117 

1619.4155 

-169.99 

-65.1936 

0.0117 

1647.9070 

-169.92 

-65.3336 

0.0118 

1875.4736 

-169.86  " 

-65.4622 

0.0118 

1902.1543 

-169.80 

-65.5849 

0.0119 

1927.9629 

-169.74 

-65.7021 

0.0119 

1952.9902 

-169.69 

-65.8140 

0.0119 

1977.2053 

-169.64 

-65.9210 

0.0120 

2000.6611 

-169.59 

-66.02  35 

0.0120 

2023.3828 

-169.54 

-66.1216 

0.0121 

2045.3960 

-169.50 

-66.2155 

0.0121 

2066.7292 

-169.46 

-66.3057 

0.0122 

2067.4021 

-169.42 

-66.3921 

0.0122 

2107.4360 

-169.38 

-66.4751 

0.0123 

2126.8579 

-169.34 

-66.5548 

DESCRIBING  FUNCTION  DATA 
FREQUENCY  =  39.177 


0.0089 

0.0090 

0.0090 

0.0091 

0.0091 

0.0092 

0.0092 

0.0093 

0. 0093 

0.0094 

0.0094 

0.0095 

0.0095 

0.0096 

0.0096 

0.0096 

0.0097 

0.0097 

0.0098 

0.0098 

0.0099 

0.0099 

0.0100 

O.OLOO 

0.0101 

0.0101 

0.0102 

0.0102 

0.0103 

0.0103 

0.0103 

0.0104 

0.0104 

0.0105 

0.0105 

0.01 06 

0.0106 

0.0107 

0.0107 

0.0108 

0.0108 

0.0109 

0.0109 

0.0110 

0.0110 

0.0110 

0.0111 

0.0111 

0.0112 

0.0112 


ALPHA  =  0.0089 


N 

1187.9441 
868.9451 
629.1567 
434.7249 
278.2578 
177.3013 
181.9892 
265.9189 
368.7612 
472.1443 
571.8198 
666.8135 
757.0452 
842.7185 
924.1157 
1001.5293 
1075.2397 
1145.5071 
1212.5645 
1276.6306 
1337.8967 
1396.5391 
1452.7180 
1506.5728 
1558.2449 
1607.8503 
1655.5020 
1701.3018 
1745.3438 


1787 
1828 
1867 
1905 
1942 
1977 
2011 
2043 
2075 


7153 
4968 
,7668 
,5940 
,0425 
.1741 
.0481 
.7129 
2268 


2105.6323 

2134.9736 

2163.2981 

2190.6414 

2217.0413 

2242.5378 

2267.1589 

2290.9424 

2313.9160 

2336.1084 

2357.5510 

2378.2678 


PHASE 


22.79 
25.80 
30.09 
37.20 
51.26 
83.21 
130.27 
156.81 
168.57 
174.70 
178.40 
-179.14 
-177.39 
-176.08 
-175.06 
-174.24 
-173.57 
-173.01 
-172.53 
-172.12 
-171.77 
-171.45 
-171.17 
-170.93 
-170.70 
-170.50 
-170.32 
-170.15 
-170.00 
-169.85 
-169.72 
-169.60 
-169.49 
-169.38 
-169.28 
-169.19 
-169.10 
-169.02 
-168.94 
-168.87 
-168,80 
-168.73 
-168.67 
-168.61 
-168.55 
-168.50 
-168.45 
-168.40 
-168.35 
-168.30 


-1/N(DB) 


-61.4959 

-58.7798 

-55.9752 

-52.7643 

-48.8889 

-44.9742 

-45.2009 

-48.4950 

-51.3349 

-53.4815 

-55.1452 

-56.4801 

-57.5824 

-58.5136 

-59.3145 

-60.0133 

-60.6301 


-61 
-61 
-62 
-62 
-62 
-63 
-63 
-63 


1800 
6741 
1213 
5285 
9010 
2436 
5593 
8527 
-64.1249 
-64.3786 
-64.6156 
-64.8376 
-65.0460 
-65.2419 
-65.4264 
-65.6006 
-65.7652 
-65.9209 
-66.06  85 
-66.2084 
-66.3413 
-66.4676 
-66.5879 
-66.7023 
-66.8114 
-66.9155 
-67.0148 
-67.1096 
-67.2003 
-67.2869 
-67.3699 
-67.4492 
-67.5252 
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DESCRIBING  FUNCTION  DATA 
FREQUENCY  =  45.000 

ALPHA  =  0.C075 


B 

N 

PHASE 

-1/N(DB) 

0.0075 

1644.6  440 

25.95 

-64.3214 

0.0076 

1229.9536 

29.06 

-61.7978 

0.0076 

918.9026 

33.25 

-59.2654 

0.0077 

667.0461 

39.63 

-56.4831 

0.0077 

462.0662 

50.64 

-53.2941 

0.0078 

311  ,9058 

71,84 

-49,88U5 

Co  0078 

253 o 9 309 

108.21 

-48  .,094  3 

0.G079 

306.6  291 

141.56 

-49.7379 

0.0079 

414.2766 

159.66 

-52.3458 

0.0080 

535.0110 

169. 16 

-54.5673 

0.0080 

655.6736 

174. 74 

-56.3337 

0.0081 

772.2756 

178.35 

-57.7554 

0.0081 

883.o467 

-179.12 

-58.925b 

0.0082 

989.5684 

-177.26 

-59.9089 

0.0082 

1090.1665 

-17  5.84 

-60  .7498 

0.00  8  3 

1165.697  3 

-174.71 

-61.47  9  5 

0.0083 

1276.4556 

-173.79 

-62.1201 

0.0084 

1362.7395 

-173.03 

-62.6882 

0.0084 

1444.8337 

-172.39 

-63. 1964 

0.0085 

1523.0073 

-171.84 

-6  3.6540 

0.0085 

1597. 5042 

-171.37 

-64.06  8  8 

0.00  86 

1666.5491 

-170.95 

-64.446  8 

0.0086 

1736.3530 

-170.59 

-64.7928 

0.0087 

1801.1094 

-170.26 

-65.1108 

0.00  87 

1862.9915 

-169.97 

-65.40  42 

0.0086 

1922. 1577 

-169.70 

-6  5.67  5  8 

0.0088 

1976.7661 

-169.46 

-65.9279 

0.0089 

2032.9443 

-169.25 

-66.1625 

0.0089 

2084.6  2  64 

-169.05 

-66.3614 

0.0090 

2134.5244 

-168.86 

-6  6.586  0 

0.0090 

2182.1536 

-168.69 

-66.7777 

0.0091 

2227.6096 

-168.54 

-66.9576 

0.0091 

2271.5955 

-168.39 

-67. 1266 

0.0091 

2313.5918 

-168.25 

-67.2857 

0.0092 

2353.6892 

-168.13 

-67.4357 

0.0092 

239  2.5  610 

-168.01 

-67.5773 

0.0093 

2429.6802 

-167.90 

-67.7110 

0.00  9  3 

2465.3181 

-167.79 

-67.6375 

0.0094 

2499.5369 

-167.69 

-67.9572 

0.0094 

2552.3938 

-167.60 

-68.0706 

0.0095 

2563.9534 

-167.51 

-68. 1782 

0.0095 

2594.2688 

-167.43 

-68.2803 

0.0096 

2623,3877 

-167.35 

-68.3772 

0.0096 

2651.3550 

-167.27 

-68.4693 

0.0097 

267  6.2  290 

-167.20 

-68.5569 

0.0097 

2  704.0447 

-167.13 

-68.6403 

0.0098 

2728.6430 

-167.06 

-68.7196 

0.0098 

2752.6655 

-167.00 

-68.7951 

0.0099 

2775.5508 

-166.94 

-68.6670 

0.0099 

2797.5315 

-166.88 

-68.9355 
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APPENDIX  B 

Computer  program  for  simulation  of  the  system  and 
numerical  data. 
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###♦* ' **#♦* ###£$. ^^J^** $$$:$$ ***^c^;. 

C(1)=MAXIMUM  VALUE  OF  INPUT  SINE  WAVE. 

C(2)=VALUE  OF  MOTOR  ARMATURE  INDUCTANCE. 

C(3)=VALUE  OF  MOTOR  ARMATURE  RESISTANCE. 

C(4)=VALUE  OF  MOMENT  OF  INERTIA. 

C(5)=REFERENCE  VOLTAGE. 

C(6)=GAIN  OF  AMPLIFIER  BEFORE  DEAD  ZONE. 

C(7)=MAGNITU0E  OF  DEAD  ZONE. 

C(8)=GAIN  OF  AMPLIFIER  AFTER  DEAD  ZONE. 

#:$:**# #  :£  #  #  £ ^c^^^:^ £#:$:;»;# ##:£;{;# #  #  #  t}c  # . 

DIMENSION  X(30) ,XDOT(30) ,C( 15) 
C(10)=l. 
IBP=1 

PI  =  3. 141592 
X( 101=0.0 
X(8)=IBP 
1  CALL  INTEG2(T,X,XD0T,C) 
X(9)=X(  D-97.0 
X(10)=T-2.9485 
XMAX  =  C(  l)*SIN(PI/3. ) 
THA=2.0*PI*60.0*T 
PHASEA=ABS(C( 1 )*SIN(THA) ) 
PHASEB=ABS(C(l)*SIN(VHA+2.0*PI/3.) ) 
PHASEC=ABS(C( 1 ) *S I N ( THA+P I / 3. J ) 
X(6)=AMAX1(PHASEA,PHASEB,PHASEC) 
IF(PHASEA.GE.XMAX)IPH=1 
IF(PHASEB.GE.XMAX) IPH-2 
IF(PHASEC.GE.XMAX) IPH=3 

206  ERROR=0.318*X(1I-C(5) 
X(8)=IPH 
VK1=C(6)*ERR0R 

IF(VK1  .GE.C(T)  )  VDZ=VK1-C(7» 

IF(VK1.LT.C(7) )  VDZ=0.0 

IF( VK1.LT.0.0)  VDZ=VK1 

VK2=VDZ*C(8) 

GO  T0(  100,101,102, 103, 104) , IBP 

100  CONTINUE 

IF(VK2)    207,207,208 

207  IBP=1 
VIN=X(6) 
X(5)=VIN 
X(7)=0. 

X(3)=(1.0/C(2) )*<VIN-0.95*X(1)-C(3)*X(2)  ) 

XD0T(2)=X<3) 

IF(X(2).LT.-0.01)     X(2)=-0.01 

X(4)=(  1.0/C(4))*(0.95*X(2)-0.086*X( 1) ) 

XDOTQ  )=X(4) 

VLAST=VIN 

X(8)=IBP 

GO    TO    1 

208  IF( VIN.LT.VLAST)     GO    TO    207 
GO    T0( 101 f 102, 103) , IPH 

101  VIN=PHASEA 
IF(VK2)     207,207,105 

105  CONTINUE 
X(7)=l. 
X(5)=VIN 
IF(VIN-. 75)300, 301, 301 

300  VIN=0.0 
X(5)=VIN 
IBP=5 

GO  TO  302 

301  IBP=2 

302  X(3)=(1.0/C(2) )*(VIN-0.95*X(1)-C(3)*X(2) ) 
IF(X(2)  .LT.-O.Oi)  X(3)=0.0 
XD0T(2)=X(3) 

IF(X(2  )  .LT.-O.OI)  X(2)=-0.01 
X(4)=(1.0/C(4) )*<0.95*X(2)-0.086*X( 1) ) 
XD0T(1)=X(4) 


7^ 


102 
106 

303 


304 
305 


103 
107 

306 


307 
308 


104 


309 

310 


VLAST= 

X(8)=I 

GO  TO 

VIN=PH 

IF( VK2 

CONTIN 

X(7)=2 

X(5)=V 

IF(VIN 

VIN=0. 

X(5)=V 

IBP  =  5 

GO  TO 

IBP  =  3 

X  (  3  )  =  ( 

IF(X(2 

XD0T(2 

IF(X(2 

X  (  4  )  =  ( 

XDOTQ 

VLAST= 

X(8)=I 

GO  TO 

VIN=PH 

IF( VK2 

CONTIN 

X(7)=3 

X(5)=V 

IF<VIN 

VIN=0. 

X(5)=V 

IBP  =  5 

GO  TO 

IBP  =  4 

X(3  )  =  ( 

IF(X(2 

XD0T(2 

IF(X(2 

X(4)=( 

XDOT(  1 

VLAST= 

X(8)=I 

GO  TO 

VIN=0. 

X(7)=4 

X(5)=V 

X(3)=( 

IF(X(2 

XD0K2 

IF(X(2 

X(4)=< 

XDOTQ 

IBP  =  5 

IF(VK2 

IBP=1 

CONTIN 

X(8)=I 

GO  TO 

END 


VIN 

BP 

1 

ASEB 

)  207,207,106 

UE 

IN 

-.75)303,304,304 

0 

IN 

305 

1.0/C(2) )*(VIN-0.95*X(1>-C(3)*X(2)  ) 

).LT--0.01)  X(3)=0.0 

)=X(3) 

J.LT.-0.01)  X(2)=-0.0l 

1.0/C(4) )*(0.9  5*X(2)-0.086*X( 1) ) 

)=X(4) 

VIN 

BP 

1 

ASEC 

)  207,207,107 

UE 

IN 

-.75)306,307,307 

0 

IN 

308 

1.0/C(2) )*(VIN-0.95*X(1)-C(3)*X(2)  ) 

J.LT.-0.01)  X(3)=0.0 

)=X(3) 

).LT.-0.01)  X(2)=-0.01 

1.0/C14) )*(0.95*X(2)-0.086*X( 1) ) 

J=X(4) 

VIN 

BP 

1 

0 

IN 

1.0/C(2) )*(VIN-0.95*X<1)-C(3)*X(2) ) 

).LT.-0.01)  X(3)=0.0 

)=X(3) 

).LT.-0.01)  X(2)=-0.01 

1.0/C(4) 1*(0.95*X( 2)-0.086*X( 1 ) ) 

)=X(4) 

)  309,309,310 

UE 
BP 
1 
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INPUT  DATA  NECORD 


ORDER    CF    tQUATl! 

DNS 

=       2 

INITIAL    TIME 

=    0.2948E    01 

FINAL    TIME 

=    0.40GC £    {  1 

STEP     SIZc 

=    0.8! C0E-C4 

THE    NGN- ZERO    CCI 

<STl 

'NTS,     C(  I),     AK 

C( 

1  ) 

=    0.115 Jc    03 

C( 

2) 

=    0.80G0E-02 

C( 

I  ) 

=0.80         -01 

C( 

':  ) 

=    0. 1360c     01 

C( 

5  ) 

=     D  .  3 1 0  0  E    0 2 

C( 

S) 

=    0.  l'JOOE    0  2 

C( 

8) 

=     0.1   ■'  -■•  0  E     !  2 

THc  NCN 


i-7r 


RO  INITIAL  CONDITIONS  ARE 
X  (  1  )  -  C .  5  7  5  0  E  (.2 
X(  2)  =  -.1000E-01 


THE  COLUMN  HEADINGS  AND  THE  CORRESPONDING  VARIABLES  ARE 


1  I  -V 
V  1 N 
VCUT 
I  AR 


X( 
X( 
x  ( 

X( 


G) 
5) 

1  ) 
2) 


THE  INDIVIDi 


APH  TITLES  ANC  THE  CORRESPONDING  VARIABLES 


OUTPUT  VS  TIME 
INPU1  VS  TI ME 


X(  9)  VS.  X (  1 C  ) 

x (  b )  vs.  x ( i  r : ) 
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HASAN  KOCaOGLU. SPEED  CONTROL  OF  A  DC . MOTOR. NO. 39 


TIME 


VI  N 


VOUT 


IA« 


0.29485E  CI 

0. 29501:  CI 

0.29517b  01 

0.  29532'  '  i 

0  .  2  9  5  4  E  <  1 

0.29564E  01 


0 .  2  -. 


1 


0.295 96 i  01 
0.29612  .1 
C  .  2  9  6  2  7  E    01 

0. 29643 E  i   1 

C  .  2  9  6  5  9  E  01 

0.296  75c  01 

0  .  2  9  6  9 1  £  1. 1 

0  .  2  9  7  C  7  -  L  1 

C.2S720,.-  01 

0/297  38  E  01 

(-.297!='.-  CI 

0.29  7  70'J  CI 

0.  2W56i  01 

0  .  2  9  8  C  2  '•  C  1 

0.29811  01 

). 296  331  01 

0.29849c  01 

0.29  ■:  t  :  -  0  1 

0.29881 :  0 1 

0.29897!  1 1 

0.2  •" 12 E  1 

0 . 2  9  9  2  i  01 

0. 29944 E  C 1 

0. 299  60c  01 

0 • 2  9 9  76c  0  1 

0  .  2  9  9  9 2  E  01 

0.30007E  01 

0  •  3  0  0  2  3  E  01 

0. 30039c  ..  1 

C.30''55::  CI 

0.300 71 E  01 

0'.  30007c;  01 

0. 30102E  01 

0 . 30 113  0  01 

0.3O134E  01 

0.30150c  01 

0.301 66 E  01 

0.30182E  CI 

0 . 30 197c  C 1 

0.30213E  01 

0.30229-  01 

0.30245E  CI 

0.30261c  01 


0.11439c  Oi 

0 .  9 7 8 3  2E  0  2 

0.1 1439 E  t3 

O.lOal^F  03 

0 . 1 1 1 44  E  .  3 

0.11269c  0  3 

C  .  1  36  LIE  0  3 

Z  •  6  2  £  -V  2  c  C  2 

0  .  0 

0.1 145 IE  0  3       0.9743 3 E    02 


0. 97497 £  0  2 

0. 97489c  02 

0.9 743 2E  0  2 

0.97478E  02 

C .     7476E  02 

0 . 974  73  E  0  2 

0  .  -  V  48  4  c  02 

0 . 9  7  4  9  2  E  02 

3 ;:  0  2 


^c 


0.10 75-0  0  3 

0 . 1 1 1 7  5  E  0  3 

0. 11243E  C 3 

0.1  j  6  j>  y  E  0  3 

0.114  75E  03 

3.99150c  0  2 

0.49255c  u2 
0 *  0 

' : !  1  i  2  i  o  r  0  3 

0 . 1 0  7  0  6  £  03 

U.11465E  0  3 

0.99797b  0  2 

0.11472E  03 

0. 106 77 E  0  3 

0.112320-  0  3 

C  .790  >  IE  02 

D .  1 S  4  2  a  E  U  2 
0.0 


D  . C 

3 . 0 

0.  10629E 

0  3 

0.  112  59-t 

0  3 

0.  1H36E 

03 

0  •  i  u  7  9  6  E 

u) 

0  .  1  1  s  >  '♦  c 

03 

0.  1C10  3E 

0  3 

0  . 1 1 4  8  7  E 

03 

0  .  9 1 9  ii  3  E 

0  2 

0. 3724  IE 

02 

0  .  0 

0 . 0 

C  .  0 

0.0 

0 . 1 1 49 2 E  03 
0.10  530E  0  3 
0.113u3E  o 3 
0  .  1 1 0  9 1 E  '.3 
0.10383c  U 3 
0.11 4 1  6  E    0  3 


0.  -7475 E  0  2 

C  .     7  4  71^  02 

0.9747  IE  02 

0  .     7 4  7  4 1  0  2 

0.97480E'  02 

0  .  9  7  4  9 1  £  02 

0  . '?  7  3  ^  2  ^  0  2 

o.r:o:.'.v  02 

0.97493E  02 

0.  574?  3!-  02 

G. 97476c  0  2 

0. 57472E  0  2 

C  .  9  7  4  7  3  c  J  2 

/ .  .  7 +76 :  o2 

•j  .  c  7  *  3  4  E  0  2 

C.97494E  02 

0  •  9 7 5  08  E  0  2 

1  .  97  5 16c  02 
0.9751? :-  0  2 
0.97504E  0  2 

0 . 97494E  02 

0.97434E  02 

0.97476E  02 

0.974 7 2E  0  2 

0.9747  1v  0  2 

0 .  37  0  7  =•;•  c  0  2 

0.9748  IE  02 

0.97491E  02 

0.^75'>-:  02 

0.97517E  02 

0.97  5  170  0  2 

0.97507E  02 

0.97493E  02 

0.97 43  HO  0  2 

0.9747 9 E  02 

0.97473c  02 

O.97071E  02 

0.97473E  02 

0.97478E  u2 

0.97486E  02 


•0.  L0000E-0.1 
0.33 633 E  01 
0.2  82  3  5E  01 
0.56674E  01 
0.96778E  01 
0.123140  02 
0.16117c  0  2 
0.14752c  02 
0.2340OE  01 
0.26  836E    00 

0.31946E  01 
0  .72116c  0  1 
0 . 9  8  8  4  2  E  01 
0.13750E  02 
0.16  537E  02 
0.19769E  02 
0.16242E  0  2 
0,  12784E  C  1 
•0. 100 00 E- 01 
0.34966E    00 

0.43853E  01 
0 . 72  9 1 7  E  01 
0.10  700  2  0  2 
0.14086*  02 
0.16  349E  02 
0.206l2t  02 
0.21466E  0  2 
0.12  6  97E  0  2 
■0. 100  00E-0] 
■0.  1Q000E-01 

•0. 100  00E-01 


-0.  10  00: 


)E- 


01 


t;.  37,,  2  30  01 

0.643  25E  01 

0.104C6E  02 

0.13  106E  0  2 

0 .  16  6  1  00  0  2 

0.197000  02 

0.22220E  02 

0 . 1 6  6  84  E  0  2 

0.67179E-01 

•0.100  0  IE- 01 

■0.100COE-C1 

■0.10000E-01 

0.21703E    0  1 

0.52198E    01 

0.90 96 8E    01 

0.11 7  510    0  2 

0.156 600    0  2 

0.182390     02 
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HASAN    KOCAOGLU.  SPEED    CONTROL    OF    A    DC. MOT  0 *.N0.39 


TIME 


V  I  N 


VCUT 


I  A  R 


o 

,  3C  277E 

CI 

0  .  1 0  2  2  2  £ 

C  3 

0.97  4  9  9  i 

0  2 

0.21747E 

0.2 

0  1 

►30292! 

01 

C  .  3  •>  9  5  4  E 

02 

0.975 12E 

0  2 

0.19085E 

02 

•J 

> 30308 

01 

0. 

0 .  V  7  3  1  ?  E 

Q2 

0.50 30 IE 

01 

c 

.  3u32<  t 

:  1 

G  .  0 

0.97! 

62 

-0.  100C0E- 

-01 

u, 

.  3  v '  3 

01 

:.  3 

0.9/  4  9  £  .. 

0  2 

-0.100(  0 E- 

-(  1 

►  303:  b  i 

t  1 

'. .  0 

0.974! 

0  2 

-0.  10000  E- 

-01 

c 

.  30372c 

01 

i  .  114J0E 

u3 

0.974 

0  2 

0. 18 24  IE 

01 

,  303  6  "^  r 

01 

0. 10 23  IE 

u3 

0.9747 3: 

02 

3. 55 8 9  IE 

01 

c 

,  3  Q  /  : 

01 

(    .11-,: 

(  3 

1  .   :  7  4  7  1 : 

02 

0.87767E 

01 

/  ■ 

►304  3     . 

01 

.;  .  i  i423E 

03 

0  .  -  7  4  7  3  E 

02 

0. 117 98E 

0  2 

Oi 

►  304 

01 

0. 1 1352E 

(  3 

0.97  ^  7 8 E 

J  9 

C. 15  507E 

0  2 

c 

,30431  . 

01 

0 . 1 1 C 226 

0  3 

0.:'?4-7:: 

0  2 

0 .  1 8  0  7  8  E 

0  2 

0. 

,  3 1      6  /  E 

Cl 

0  . 1 0  9  6  2  c 

0  3 

0.  ?74<  =  E 

9  2 

0 . 2 1 8  94  E 

0  2 

0 

:  2 

i  1 

0.7 1112c 

0 .  3  7514  E 

02 

0.21640E 

0  2 

0, 

,  30^9H 

..1 

1    .  8041 3 t 

ul 

0.97523E' 

02 

0. 10  9  95E 

0  2 

0 . 

,305]      : 

01 

o.o 

0.97517E 

02 

-0.100      il 

-01 

c  . 

,"3C  5 : 

'   1 

0  •  u 

0 .  9  7  5  0  7  £ 

02 

-0.  10  0  00  E- 

-01 

305     ■ 

;  1 

-  . 

0. 97497 E 

0  2 

-0. IOC  00  2- 

-01 

i  ( 

01 

0 . 974  8  7  E 

02 

-0. 100 00 E- 

-01 

0 

►  3051 7 

1 

0 .  1 1 0  00  E 

"3 

... 374 78E 

0  2 

0.  24-    ■ 

01 

a  (   »~  c  3  r 

01 

0.  I  136  .■'- 

•j  3 

.974  7  ';  ': 

0  2 

0.92339.- 

01 

0  i 

09  E 

0] 

C  .  1039  2c 

0  3 

0 . 974  7 1 E 

02 

0.9 0 2 8 3 ! 

01 

sJ   i 

, 

01 

O.liSOoE 

0  3 

0  . 9  7  4  7  2  E 

•'2 

0.12 0 7 7 : 

0  2 

.  . 

,  30c-.  ■   ]  , 

01 

0 .  1 0  3  I  •• 

C  3 

.  .     74 78 E 

02 

0 .151;    E 

j  2 

^.  . 

, 30656c 

01 

0.113: 

03 

0.  974 

0  2 

0.1870-9 

02 

c  < 

,  30  67i  E 

01 

■   .85  2  7-ti. 

0  2 

3.97  4 

0  2 

0.20 3 97E 

0  2 

1  , 

.  3  .  6  b  I  c 

1  1 

1,  27  J  •.  St 

C  2 

0 ,  375  £  ii  E 

)2 

0.  13134E 

02 

c 

,  3(  704  :. 

c  1 

0  .  ( 

C  .  37  5o ■ 

0  2 

-C.  100G0E- 

-  0  1 

V    « 

, 30720E 

v.    1 

0 . 3 

0  .  9  7  4  c  4  E 

0  2 

-0.100 

-  0  1 

0 

,30736c 

01 

0 . 0 

0.974 

02 

-0.  100  C(     - 

-',1 

0  < 

,  3  0  7 :  1 E 

31 

0.1 0 2  5  6  E 

0  3 

0  .9  7  4  7"/  E 

02 

0.  301  »0E 

0  1 

u. 

, 307  67; 

01 

0. 11^07E 

u3 

0 .  9 74 7 2  : 

0  -> 

0  ■  6  7  4  8  3  E 

01 

0. 

1  3  "J  7  &  'i  i 

01 

0  .  10  30  or: 

..3 

0  .  9  7  4  7  1 E 

0  2 

0.950*  '■'■: 

01 

Oi 

, 307 99E 

01 

0.11 u  1 1 1 

j  3 

0 .  9  "7  4  7  ->  E 

i    ? 

0.13456E 

02 

0, 

,  3  C   J 1 5  E 

01 

0 .11321  & 

03 

C. 97480= 

0  2 

C.16048E 

0  2 

u  < 

,30? 3 It 

Ul 

U  .  1  0  -+  "7  v  E 

03 

0.9749  IE 

0  2 

0. 3  97  30c 

9  2 

.';. 

, 3CS46E 

01 

0.60^8 3 1 

02 

0.975 02 E 

0  2 

0.17950E 

0  2 

0, 

, 30  862c 

c  1 

(  .0 

0 .  9  7  5 0 6 E 

02 

0.49 03  IE 

01 

0, 

303  7  E E 

01 

0.  3 

0.97497 1 

02 

-0 .  10  3  0GE- 

-  0  1 

Oi 

308 94 E 

ul 

0 .  0 

0.  0-74  5  7'c 

02 

-0.  10  000  E- 

-01 

0, 

,309 10 E 

01 

0  .  1  1  1 J  5  E 

0  3 

G.9747SE 

02 

0.29601E 

01 

^  « 

3C  926? 

Jl 

0. 11298E 

03 

U.97473E 

02 

0.57  1  17 E 

01 

0  < 

, 3C 94  IE 

Cl 

0 .  10550 E 

03 

0.97472E 

02 

0 .  9  5  8  3  9  E 

01 

V,   1 

309  57c 

01 

0.11  -.-  30  E 

0  3 

0.974  74E 

0  2 

0 . 1 2  5 1 5  E 

0  2 

o, 

309 73E 

L  1 

0. 1C138E 

03 

0. 97480E 

02 

0. 15 j99E 

02 

I 

, 30 9 89 E 

Cl 

C . I  1437 E 

0  3 

j. 97489 E 

02 

0. 1914UE 

0  2 

r-. 

v.     1 

310051 

01 

v.1  •  87  '.   i  -  i 

02 

0  .  9  7  5  C  2  E 

0  2 

0.21127c 

02 

u . 

,  3 1 0  2 1 E 

Cl 

0  .  3  0  9  0  3  E 

0  2 

0.97512E 

0  2 

0  .  1 4  5  3  3  E 

02 

0 . 

, 310 36 E 

01 

l.O 

0.97  5 10 E 

0  2 

-0 . 10 3  0UE- 

-01 

0. 

,  31032c 

V.     1 

0.975 00 c 

j2 

-  0 .  1 0  0  00  E  - 

-01 
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HASAN  KOCAOGLU. SPEED  CONTROL  OF  A  DC . MOTOR. NO. 39 


TIME 


VIN 


VOUT 


IAR 


0  .  3 1 0  6  8  E  0  1 

0 . 2 1 0 84 E  CI 

o .  3 1  i  oo  c  :•  i 

0 .  3 1 1 1 6  E  1 

0.3113  1  .:  i   1 

0.31147E  CI 

G.31163E  01 

0.  31179c  (  1 

0  .  3 1 1 9  5  E  -  1 

0.31211E  01 

0 .31226-  CI 

0.31242c  ' i 

G.312  5EE  01 

0.31 2 74E  CI 

0.312 90 E  01 

0.31306!  CI 

(  .31301::  0 1 

0  .  3  1  ?  3  7  -  '1 

G  .  3  i  3  5  3 ;  :.  1 

0.31369;  01 

0.  31 3 85 E  01 

0.31 4 C It  CI 

0.31  '■.  16c  C  1 

0.  314 32 E  01 

0 . 3 1 4  4  &  E  CI 

C  .  3  1  .  i  ■•-.  r  01 

0.3148Ct  CI 

0. 3  In  96 E  C3 

0.31511E  v  1 

0.31C27E  01 

0 . 3 1 5  4  3  E  01 

C . 3 1 5  5  9  E  CI 

0  .  3 1 5  7  5  E  CI 

C  .  3  1 5  ?  1  E  C  1 

0 • 3 1 5  06  E  01 

0.316  22E  Cl 

0.316 36 E  C 1 

C. 316  54c  t 1 

0.316 70 E  01 

0-31686E  01 

0.  31701 E  01 

0.317 17 E  01 

0.31733 1-  C  1 

0  .  3  1  7  4  9  E  C  1 

0 . 3 1 7  6  5  E  01 

0 .  3 1 7  8  0  E  C  1 

0.  317  9c  F  '.-1 

0.31812c  01 

0 . 3 1 S 2  8  E  01 

0.31844E  01 


0  .0 

0.  1  J073E  0  3 

0.11449c  o1 3 

J.  107 7  ■  ■-  i  3 

0  •  1 1  i  6  9  t  .  3 

G.112<*6E  0  3 

0  . 1 0  6  5  <  =  0  3 

0.11  ';  7  / '  03 

C. 990 27 I  0  2 

0.49036E  0  2 

0.17837E  0  2 


0.97490E  02 

0. 97431 E  0 2 

0.  ?7473E  02 

0.974  7  J E  0 2 

0.     7vO'.i  02 

0.'  7-   7  :•  :  0  2 

C .-  74 ?0E  02 

0.97491  -  02 

0.975 05 E  0  2 

C  .  ?7520c  02 

0  2 


0  •  7  6  c  , 

0  2 

0 .  ?  7  5 1 6  E 

0  2 

j. 1122  It 

0  3 

C .     7  5061 

0  2 

0.  10697E 

03 

0 .     f  4  9  6  E 

0  2 

0.64741c 

0  2 

0. ?74 J6E" 

02 

0.996&2E 

02 

0  .  9  7  4  7  S  E 

02 

0.11470c 

0  3 

0 .  -v  7  <i-  7  2  E 

02 

0.106? 

C3 

0 .97471 E 

02 

O.il?, 

j  3 

0.97473E 

32 

> .  1 1 1  i  2 

C  3 

0.  -)7479c 

02 

0  •  1 0  7  'i 

...3 

0.9748* 

02 

0.6  578  ii 

0  2 

0.975      )E 

02 

0. 1410 

(  i 

c .  v  7  5 

0  2 

o.O 

0  .  9  7  4 

J  2 

0 .  0 

' . 57466 E 

02 

0  .  1 1 2  5  3  E 

0  3 

0.  97477c 

0  2 

0  •  1 1 1 6  2  £ 

0  3 

0 .  ?  7  4  7  2  E 

0  2 

0.  1073     . 

0  3 

0.97 4 7  ]  E 

0  2 

0.11 4    • 

0  J 

0 .  9  7  4  7  3  E 

02 

C • 1 0 091c 

0  3 

0.97479E 

02 

0.1 14t    - 

03 

0. 9748 8E 

0  2 

./  .    '  L   1  O  j  C 

C  2 

0.975   lit 

02 

0  •  3  o  9  5  9"t 

02 

0.9751  '1 

02 

0  . 0 

0.97512E 

0  2 

0 .  0 

0. 97503 E 

0  2 

1.0 

0.97493E 

0  2 

0 . 1 0  1 5 1  E 

0  3 

0  .  9  7  4  3  3  E 

0  2 

0 . 11491c 

03 

0.97<<  75E 

02 

0.10539c 

0  3 

0.9747  IE 

02 

0.11304E 

0  3 

0 . 974  7  0  E 

02 

0. 11098E 

C3 

0.974 73 E 

02 

0.10  373E 

03 

0.^74^  ic 

02 

C .  1 1 4 1  9  E 

0  3 

0. 37490c 

02 

0.1C211E 

03 

0.9750 ',  E 

02 

C  .547  1 1 E 

0  2 

0  .  97 5  1 8  E 

0  2 

0.0 

C  .  9  7  5  2  4  E 

0  2 

0  .  c 

0.975 15E 

02 

0.9750c E 

0  2 

0.  J 

0  .  9  7  4  9  6  E 

02 

0  .  ) 

0 . 9  74  36 E 

02 

0,  10000E-0.1 

0.514  90E  VO 

C.41493E  0] 

0.70121E  01 

U.10072E  02 

0.  13  534E  02 

0.17 3 36E  0  2 

0.20119E  02 

0.23 2 SOE  0  2 

0.19660E  0  2 

0.507  16E    01 

•0.  1C0  0OE-O1 

0.  10000E— 01 

•0.  lOOOOE-Ol 

•0. 10000E-01 

0.25 3  0  IE     .1 

0.60  2S5E    01 

0.89137E    01 

0.12 3 09E    0  2 

0 .  1 5  3  9  2  E    ' :  2 

0.192  02c  0  2 
'.'.13224c  02 
0.64137E  01 
•0.  10  0  00  £-01 
•''■.  1000  »£-01 
0.2  93S8E  01 
0 . 5  7  2  7  5  E  01 
0. 97 0  91 E  01 
0.12465E  02 
'"  .  15  934E    0  2 

0.19 11  IE    '.2 

0.21556E    02 

0.15  983E    02 

•0.  10000E-01 

■0.100  00  E-01 

•0. lOOOOE-Ol 

0.10 0  7  IE    00 

0.3400  IE    0  1 

0.65089E    01 

0.10  3  7  IE    0  2 

0.130C5E    02 

0.168  9  2E    ,:2 

0.19  5 0  5 E    0 2 

0.22934E     02 

0.20  219E    0  2 

0.610 70 E    01 

•0. 100CGE-01 

■O.K    IOOE-01 

■0. 10000  c-01 

•0. 10000E-01 
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H,oAN  KGCAQGLU. SPEED  CONTROL  OF  A  DC . MOT  OR. No. 39 


TIME 


VI  N 


VGUT 


IAR 


0  . 

318t"  '■" 

01 

0  .  10  2  i  u  l 

0  3 

0. 

97478c 

02 

0  *  27  8  96E 

01 

<> 

31c   " 

01 

.  1  1  •'-..'  t 

03 

\J   • 

-  f472t 

jZ 

J  .  6  0  3  1  5  E 

3  1 

0. 

31  J91t 

C  1 

■0.1043: E 

03 

c. 

rTMlc 

02 

0 . 9087 0 3 

01 

0. 

3  1 ;.    u 

01 

0.113  • 

03 

97473 7 

0  2 

0. 12  846E 

'.•2 

0. 

3 1  ?  ^  3 ; 

01 

0.11023E 

03 

•-  • 

97479  ■ 

0  2 

3. 15v53£ 

02 

0. 

?!  •. 

'.  ] 

C.  1   ;  ;  -<  ?l 

03 

0  • 

9743  8c 

3  2 

0.1931 

0  2 

1. 

319     : 

CI 

0.70921E 

02 

0 . 

97  3 

02 

0.190-iih 

0  2 

V.'  • 

31 

I  1 

3.77: SSfc 

01 

0. 

975!  7  - 

)2 

3.  34  3  37  3 

01 

0  . 

3] 

01 

i .   •  3  0  '-*■  3 1 

0  2 

u  • 

97  \ 

0  2 

- 0 . 100  .  1 E  - 

-01 

0. 

32..  . 

'  i 

1 .  1  0  3  i 

j  3 

0  . 

r:    (  4  3 

02 

-0.1000  0 ;:  - 

-(  1 

0. 

32i  1  r E 

-  1 

f.-  .  1  l  3  c  3  c 

0  3 

974 

0  2 

0. 207 3 JE 

01 

0 . 

32G  3 

.  1 

0. 1099  IE 

03 

0 . 

974" 

3  2 

0.4  86  78E 

01 

1  . 

LI 

0 .  1 0  9  9  2  E 

03 

L  « 

974 72E 

02 

0.  83  :; 

(-1 

320<  . 

01 

. 1136  7E 

0  3 

0 . 

9  /  4  7  3  E 

-.2 

0 .  1 1  5  8  Ql 

0  2 

0  • 

3  2  0  6 IE 

C  1 

0  .103    21 

C  5 

u . 

37476  E' 

02 

0  .  1 5  7  o  3  3 

02 

u  . 

323c?t 

01 

J  «   Li  ,  J  i>  . 

13 

f 

97486E 

o2 

0. 162  24E 

02 

(  . 

3211: 

C  1 

u.:o^3   E 

3  2 

^  • 

97499 E 

02 

0.21O77E 

02 

•;. 

3212-j 

01 

0 .  42 9  5 <* E 

02 

u. 

97  5  1 0E 

02 

0.16 5  09  3 

0  2 

0. 

321  4£i 

01 

0  •  0 

n 

■  /   • 

r  7  5  1  1 .-. 

02 

0.62462E 

00 

(j  . 

32160i 

01 

C  .  J 

..  :  c 

97 5 v  It 

0  2 

-0.  10000E- 

-01 

0  . 

3217c  :-. 

C  1 

.  n 

0. 

974   i  1 

02 

-C  .  ICOCOc- 

-01 

'•J  • 

32192 

C  1 

0.10436c 

0  3 

0 . 

:'  7  3  -  2  E 

oz 

0.6698  2E 

00 

0  • 

3  2  2  I  = 

}1 

.  1  i'i  •/.: 

03 

0. 

-  ?  i  7  v : 

02 

0 .  3  S  3  7  3  E 

)1 

0. 

322  24c 

0} 

0  . 1 0  2  6  7  E 

I  3 

'  . 

i  i  70E 

0  2 

0. 70550c 

31 

c . 

322'.    ■: 

CI 

3  . 1 1  ; 

..3 

r 

« ■  • 

974701 

o; 

0. 10734c 

)2 

Li  • 

.       5  3  : 

I  I 

3.10  •.' )    i 

03 

L  « 

974  ' 

02 

0.  13-. 

02 

• 

3227i. 

01 

0  •  11 0 3  c  E 

0  3 

0. 

974-  IE 

0  2 

0.  17 3  133 

02 

»• 

322  87 

01 

J • 7  3  90  5 E 

L-2 

\_  * 

v7-V-  IE 

0  2 

i- .  1  7  5  3  :  - 

0  2 

32  3, 

CI 

0.11599c 

02 

c. 

9  74  <  -j  3 

02 

0.76?  03  - 

01 

0  • 

3231c  ■ 

01 

0.0 

0  • 

9743 6E 

0  2 

-0.  100303- 

-  0  1 

'J. 

3  2  3  ?L  • 

CI 

0.1021 0 £ 

0. 

97479 E 

02 

0 .  1 7 1  3 1 E 

01 

0. 

3  2  ^  : 

01 

0.11 *19E 

0  3 

0. 

974 7 2E 

0  2 

0.543  3  IE 

01 

« 

32  36£  t 

01 

0.  10S7  3-E 

0  3 

v'     • 

3747  )E 

U2 

0.32  2  52E 

01 

p. 

32362 £ 

01 

0.11  C v  c t 

0  3 

>J  ♦ 

97472E 

02 

0. 12 191 E 

02 

0  • 

323'. 

01 

0. 1130  3c 

0  3 

0. 

374 77 E 

02 

0. 147  973 

0  2 

D. 

324 14 t 

01 

0. ]  yt  .■■, t 

0  3 

0 . 

9  7  •: -3  3  E 

02 

0  .  1 8  5  2  3  L 

'.:Z 

324 3C  '. 

01 

0  .  o  1  3  -i  J  c 

02 

*  0. 

9  74'"*  3  c 

3  2 

0. 16395E 

02 

0 . 

32445E 

01 

0.0 

0. 

97499 E 

0  2 

0.41 2 7 93 

01 

0. 

3246] 

01 

'  •  u 

0. 

<=  7  's  39  : 

02 

-0.10000c- 

-01 

324 77  £ 

31 

0.  10830c 

03 

r 

9748 0 E 

OZ 

0.98985E 

(  0 

32493E 

0  1 

0.11132E 

0  3 

c « 

97473E 

k)Z 

0.50 5 92E 

01 

0. 

32 5 09 £ 

01 

J. 11279E 

0  3 

J  • 

374706 

02 

0.777233 

01 

32525 • 

01 

0  .  1 J  5  9    E 

03 

•> 

3747  LE 

02 

0.116393 

0  2 

c . 

"2  ~>  E  .',  ."    ,. 
-i  C-  .-    "".  • 

Ll 

0  •  1 1 4  8  o  E 

03 

o . 

:  7  s-  7  5  E 

02 

0. 14505c 

0  2 

0. 

3255$ c 

01 

0 • 1 00 ?OE 

03 

v    * 

57434E 

02 

0. 176 98t 

0  2 

0. 

32572 

01 

0.11  4  .  J  i 

0  3 

(.:  . 

9  7  4  9  5  E 

i.3 

C.Elv  -,1: 

0  2 

0. 

32b£ . . 

01 

i.;   •    3  8    -r 

u2 

v.   . 

375] 

0  2 

3.  231 

02 

0. 

3  2  o  0  4  E 

01 

0.31884E 

02 

3 . 

97  5  2 3 l 

02 

0.  16  6  45  E 

0  2 

0  . 

32620: 

C  1 

0. 

97522E 

•3  2 

-0.  L0QC0E- 

-01 

0. 

32635E 

01 

0.0 

0. 

9  7  3  1 3  E 

02 

-0.  10000c- 

-01 
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HASAN  KOCAOGLU. SPEED  CONTROL  OF  A  DC . MOTOR. NO. 39 


TIME 


V I  N 


VOUT 


IAR 


0.32t  ;  1.. 

01 

0  .  0 

0.97503E 

02 

-0.  10000E- 

-01 

0.  32-6  67! 

01 

l  .0 

0. 97493E 

0  2 

-0.  10000E- 

-0  1 

C.  32; 

C  1 

0.11458 ': 

03 

0.974835 

0  2 

0.  3  95  99E 

00 

0  •  3  2  6  9  9  E 

01 

0.  10  <  '>3fc 

0  3 

0  .97476b 

02 

0.33376E 

01 

0.  3271:".. 

1 1 

i . U193E 

03 

0.97472E 

02 

0.73430E 

0  1 

0.32  7. 

01 

0.112, 

C  i 

C .  9  7  4  7 1 E 

02 

0. 100 12E 

0  2 

0. 32746= 

"1 

0 . 1 06  6 

3 

0.97475E 

02 

0  •  1 3  8  8  8  E 

0  2 

''•.327:': 

:  1 

o .  1 1  s  r Oc 

03 

J. 9748 2E 

02 

0.1665 6 E 

02 

0. 32776- 

01 

0  .  -■  9  5  -<  i  i! 

J2 

i          .  "V         ?  c 

0  2 

0. 19920E 

0  2 

0.32  7  '  '.  = 

C  i 

■J    .    '.    i    7  .     . 

02 

0.9  75C3E 

)2 

0  .  1 6  5  0  2  E 

"2 

0  .  3  2  8  1 C  " 

01 

0  . 0 

0. 975  04c 

02 

0. 16445E 

01 

0.32    2 

01 

C .  0 

0  . '  7  ■',- 

)2 

-0. 10  0  00  E- 

-01 

0.3284  i 

01 

0.,' 

0.97465E 

02 

—  0  •  1 0  0  0  0  E  - 

-01 

0  .  3  2 e '.  (  c 

01 

j  .  10734c 

0  3 

0  .  9 7 4  77  E. 

0  2 

0.3  7  334E 

01 

0. 328  73s 

f  1 

0.11  -';■  ■ 

0  .  9  7  4  7  3  E 

02 

0  .  cO'JiZt 

0] 

(  .32889 

"1 

0  .  1 0  0 1 7  E 

C3 

0.    '*/  ',  722 

02 

0. 10036c 

0  2 

(.: .  329 J    . 

01 

0.1147 7  b 

■'3 

C  .9  74  7 5  E 

02 

0.134 19c 

0  2 

! . .  ?.  ?  -"  ■  2  >  i 

0] 

0.  10  6  *  3E 

03 

.-»      OT£i  '•'■  7  !- 

0  2 

{^.1621'. 

0  2 

0.3293 • 

01 

0.112 

0  3 

0.9  7  ■■-':■  t  2  E 

0  2 

0. 19970E 

0  2 

o .  3  2  9  1  2  : 

01 

U  .  ;  v  z  J  i  '; 

02 

;  .  vv;-  »5E 

02 

0  .  2  0  9 1  o  £ 

■J  2 

0.32  9 6 i  E 

01 

0.  19  1865 

C .  9  7  5 1 4  E 

02 

0.  12236E 

0  2 

3.  329*  ■  E 

01 

0  .  I  ! 

0  .  'i  i  :- 

0  2 

-0.  IOC  0  >E- 

-01 

C-.3; 

£  1 

C  •  J 

02 

-0.10CQG 

-01 

0  .  3  3  0 1 1  E 

0  1 

( ■ .  c 

0.  974 90£ 

0  2 

-0.10, 

-  0  1 

J  •  3  3  L  3 1  E 

C  1 

0 . 1 0600c 

0  3 

0.97430c 

0  2 

0.61 872 E 

00 

81 


INPUT  DATA  RECORD 


ORDER  Dc  EGUATIGNS  =  2 

INITIAL  TIME  =  0.2 59 ^E  01 

FINAL  TIMt  =  3. 3500 E  01 

STEP  SIZE  =  0.800UC-04 

THE  NGN-ZERO  CONSTANTS,  C(I),  &RE 


TH 


MON-ZERi 


C( 

1) 

= 

Q.115JE  J3 

C( 

2) 

= 

!.- .  SO yju E - 0 2 

C( 

3  ) 

= 

..•>...'  -CI 

C( 

<  ) 

= 

0.138UE  01 

C( 

5  ) 

= 

0  •  3  2  8  u  c  3  2 

C( 

) 

= 

;..' .  I  3  0  J  E  0  2 

C( 

i  ) 

= 

C.10C0E  02 

IN 

I  T  I  A  L 

C  0 K  C I T 1 0 ;  N S  ARE 

X( 

1  ) 

= 

,  .  L  0  3 1 E  0  3 

X( 

2) 

= 

0  •  14  2  -<  c  0  2 

THE  COLUMN  HEADINGS  AND  THE  CORRESPONDING  VARIABLES  ARE 


TI^E 
V  I N 
VCU7 

i  AR 


X( 

0) 

X( 

5) 

X( 

I  ) 

X( 

2) 

TH, 


INDIVIDUAL  GRAPH  TITLES  AND  THE  CORRESPONDING  VARIABLES 


0 IT  PUT  VS  TIME 


INPUT  VS 


I  " ': 


X (  c  )  VS.  X  <  1 0  ) 

X (  b )      VS.  X (  1  J ) 
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HASAN  KOCAGGLU. SPEED  CONTROL  OF  A  DC . MOTOR. NO. 59 


TIME 


VI  N 


VOUT 


IAR 


,  ?^3TH 

01 

G .10927 E 

03 

;    • 

L  0  3 1  3  E 

03 

0.142  4  2E 

0  2 

1. 

, 25  955E 

0  3 

0 • 113  98c 

03 

0. 

L03  3  *E 

0  3 

0. 15  309E 

0  2 

0 

,25971c 

:  i 

0.102 S 7 E 

0  3 

w    * 

L  0  3 1 4  E 

03 

0. 18289E 

02 

c 

,259861 

01 

0  .  5  6 1 8  1 E 

02 

w    • 

!.  0  315c 

0  3 

0.14810b 

02 

0. 

, 260 02 E 

£  1 

0.0 

0. 

] 0 315  6 

03 

-0.99193E- 

-02 

u 

,  2  6  0  1 8  E 

1 

0  .  1 13^4-E 

03 

0  . 

1  0  3  1  -  •: 

0  3 

0.410  07E 

0  0 

c , 

,  2  6  C  3  4  E 

01 

1.11017c 

03 

0  . 

L  0  3 1  3  E 

03 

0.21-3696 

01 

J 

,260!      . 

,  1 

.  i 

c « 

!      313c 

03 

0.31631c 

01 

0 

,26066 E 

C  1 

'    .  11380E 

• 

L0312 E 

03 

0.68  584E 

0 1 

'J, 

,26C 6  IE 

C  1 

0. 10344E 

03 

0. 

L0312E 

0  3 

0.95  243E 

01 

j 

,2609? i 

01 

0.1150  >c 

(  j 

c . 

1 531 2E 

0  3 

0.  11536c 

0  2 

v.     . 

,261  1  3 

01 

u  .  1  -J  3  d  3  £ 

0  3 

u  * 

'  3 1  3  E 

03 

0.1 349  IE 

0  2 

0 

.26129E 

(.  1 

0.  113  7 3c 

03 

c. 

1 0  3 1  3  E 

0  3 

0.  16 06 96 

0  2 

0 

.26  Kit 

01 

C  .  1  j-/  79 E 

0  3 

s.'     • 

L0314-E 

0  3 

0.  1757^0 

0  2 

I) 

,2616] : 

01 

0. 110( 

03 

0.  ] 

L  0  3  1 5  E 

u3 

0.20  33SE 

0  2 

0 

,2617c 

f  1 

'  .72251 E 

12 

c. 

L0  316E 

03 

0.  192 02E 

0  2 

(.  , 

, 261  ;2E 

01 

0. 94651 E 

01 

0. 

L  0  3 1 7 1 

0  3 

0.77  93  IE 

01 

3. 

, 26  2CEE 

01 

0  .  5  6  5 7  Ot 

02 

c .. 

L0316E 

03 

-o.  ioocoe- 

-01 

0 , 

,262  2'  E 

01 

C  .  10  307c 

,    j 

*  J  * 

L03  3  >  E 

03 

-0. 10000 E - 

-01 

c . 

26240 1: 

0  i 

1   .  I  l  3  ',  2  c 

0  3 

'■*.    ft 

LO  3 1  4  E 

C  2 

0.16+27E 

01 

0, 

,26256c 

01 

J  .    1  w   ?  4i  c 

f '    -*■ 

0. 

L0313E 

03 

0.35-576" 

01 

262? j  E 

f.  1 

1   .Hi 

03 

\j  • 

_  3  j  3  .: 

0  3 

0.65 7 78E 

01 

2  6  2  8  7  E 

I  ] 

0.1134] 

0  3 

* 

. :  312c 

0  3 

0  .  6  2  2  4  5  E 

)1 

c . 

,263 0 ?  E 

C  1 

0  .  10  *  5  2: 

C  3 

V      • 

l  0  312? 

0  3 

0.  10^39E 

6  2 

1 

26319E 

01 

C  .  1  ]    ,  ~  v  c 

03 

\J    ft 

L0  313E 

03 

0.  1  2  ^  '■>  7  E 

0  2 

. 

,  2  6  3  3  £  £ 

t  ] 

:.  ,10250c 

l  3 

l    ft   , 

.0313c 

0  3 

0.  K87C  E 

0  2 

0, 

,26351c 

01 

0  •  1 14( 

i  ^ 

V.'    -' 

^  ft 

L  0  3 1 '  E 

03 

0. 17344E 

02 

-  « 

2  o3  6  c  c 

C  L 

( i»8  6  227c 

C  2 

L   • 

0  315? 

03 

0,  181 34E 

02 

0 

,26382c 

I  1 

0.2  6  54  5c 

■J  2 

L0315t 

0  3 

0.  lti-vt: 

0  2 

0. 

.263 96  £ 

01 

I   .0 

0. 

L0315E 

0  3 

-0. 100 OOE- 

-01 

■■•;  , 

2641  4 E 

C  1 

0.10 5 C 5 : 

03 

0  • 

103 14E 

0  3 

0.  18  060  6 

01 

,26*i  31  c 

01 

0.  1  l-',  -.  . 

C  3 

c . 

L 0 313  F 

03 

0.36182E 

01 

0, 

26446c 

01 

0.10192t 

03 

c .  . 

L 0 3  13  E 

03 

0.603  3  IE 

01 

o , 

.  26461E 

01 

C .  11424c 

03 

c. 

L  0  3 1  2 : 

03 

0. 86041E 

01 

0  . 

2  6  4  7  7  E 

01 

0  a   lu5  6  0  t 

03 

u  * 

L0312E 

0  3 

0.  10  286c 

02 

0 . 

,264  93c 

01 

U  .  1 1 1  \    -: . 

03 

0  3  1  3  E 

03 

0  . 1 3  1 5  2  E 

6  2 

0  , 

,  2  6  5 

C  1 

0  .  1 1 2<  .  - 

U  .3 

c . : 

1031 3E 

0  3 

0. 14075c 

0  2 

0. 

2  6  5  2  5  E 

01 

C  .  10557! 

03 

L0  314E 

0  3 

0. 17  346E 

02 

2  6  5  4 1  E 

C  1 

0  .11^ 5 0 z 

0  3 

0.  ] 

.031 5 E 

03 

0.19076E 

02 

1 , 

>  26556 c 

01 

0  •  9  7  7  6  3  E 

0  2 

c . 

L  0  3 1 6  E 

03  . 

0  •  2 1 0  7  3  E 

0  2 

c . 

26  572E 

01 

0  .  *  6  6 :>-  2  c 

1*2 

0.  ] 

L  0  3 1 7 1 

0  3 

0.16069E 

0  2 

c . 

.26538E 

01 

'   .0 

0 . 

L  0  3  1 7  E 

0  3 

-0. 100  006- 

-01 

n 

26604c 

CI 

0 . 0 

0  . 

L  0  3 1 6  E 

03 

-0.  10000E- 

-01 

o ! 

,2662c  E 

L  1 

C  .  0 

0.] 

L0315E 

03 

-o.  lcecoE- 

-01 

0 . 

, 26c  3tc 

C  1 

0  .  i  0  6  ( 

03 

0. 

L0314E 

0  3 

0.136626 

01 

0. 

2  6  6  5 1 E 

01 

I   . 1 1484 [ 

C  3 

0  . 

10  313  E 

0  3 

0.33369c 

01 

o, 

,  2  6  6  b  7  E 

I  1 

'■  .  3  0070E 

C3 

0. 1 

L 03 13c 

0  3 

0. 56  5  66  6 

01 

c . 

, 26683  £ 

01 

0.11450E 

0  3 

L  . 

L  0  3 1 2  E 

03 

'0.  01412c 

01 

0, 

266 99 E 

ll 

0.1O774E 

03 

0  .  . 

L0312E 

0  3 

0.98779E 

01 

, 26715E 

el 

0.  11172c 

C  3 

0. 

L0312E 

C3 

0.12725E 

02 
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rUSAN    KOCAQGLU. SPEED    CONTROL    OF    A    DC . MOTOR. NO. 59 


TIME 


VI  N 


V'jUT 


IAR 


0.26731E  : 1 

0  .  2 6  746c  CI 

0. 26762  £  01 
0. 267 7t 

0.2 67 94 J  .1 

Q.2631Cc  31 

0 .  2  6  8  2  6  E  •  I 

0.268*1:  .1 

0.2  01 

G .  2 L  6 7 3 £  CI 

0.2688'  ;  ..] 

0.269  Obi:  I 

0  .  2  6  9  2 1  £  ..•  1 

0.  2:  \  3oc  .,1 

0.269525:  i  i 

0.26968c  CI 

0.269  84-E  CI 

0.  2700<  -  (  1 

0.2701c  CI 

['• .  2  7i  '-  i  •  .  1 

o.  27^-7'  01 

u.2  7fc-*  01 

...27C7     ..  01 

0 . 2  7  0  9  5  J  01 

0.27110c  01 

(  .  27  12'.  j  Ul 

0.  27  1'  21  I  1 

0 . 2  7 1 5 c i  CI 

C  .  27174£  01 

0. 27190 E  01 

0.27205c  01 

0.2  i  c 2  1  E  C 1 

0.2723 V  E  C 1 

0.272  53E  01 

0  .  2  7  2  £  c  E  U  1 

0.27  2  6  5'i  CI 

0.273*  f  c  CI 

0.27316b  01 

0.273322  01 

0 .27348t  01 

0.273  64S  01 

0. 27380c  il 

0.27395c  11 

0.27411c  01 

0.27427:  (  1 

0.27443E  CI 

0.27-;-59E  01 

0.274 75 E  01 

0 .  2 7 4 9 0 E  CI 

C  .  2  7  5  C  c  E  0 1 


0. 11245 E  Oo 

C . 1065  5c  03 

C  .11'  1 1 E  C  3 

0.99 097E  0 2 

0.491616  02 
0  .0 
C  .  0 

I-  «  'J 

0 .10702c  C  3 

G.11467E  03 

0.9 973  IE  0  2 

.11471 :  U 3 

0.  1068  If  03 

0.11230c  0  3 

0  .  1 1 1 3  3 1  03 

0 . 1 0  74  at  03 

0.1145  7c  0  3 

0.10037c  03 

0.51481c  G 2 
0 .  •.. 

r' .  C 

0. U 793c  0  3 

0.11 4  '■*'->  E  ».  3 

0.1009:'   :  0  3 

0.1 140c"  C 3 

0.10  5  B  '*  c  0  3 

0.11282r  03 

0.11127c  0  3 

0  •  1 0  8  3  6  E  L  3 

0.11432c  03 

0.10156E  0  3 

0.53726€  02 

C  .  0 

0.111  34E  03 

0.108 78 t  03 

0.114171  03 

0.1O219E  03 

0.11496E  03 

C  .  10482c  03 

0.113291  l»  3 

0  .IK  i  uF  C  3 

0.10 S  19c  l 3 

0  •  7  '• '  0  3  3  J  ^  2 

, ,6d772E  01 

VJ  .  0 

C  .  10428  E  03 

0 . 1 1 3  5  0  E  03 

0.110  25E  03 


0 . 

1   (313E 

0  3 

0.14243E 

0  2 

L  . 

L031 VI 

0  3 

0  .  1  6  9  7  1 E 

02 

0 . 

r   314E 

03 

0. 18  642E 

0  2 

0  . 

t 0316 E 

0  3 

5.20768c 

02 

c. . 

1 0  3 1  7 1 

0  3 

C. 16143E 

0  2 

0 . 

1031  7  E 

0  3 

0.98521c- 

-01 

.  0  3  1 6  E 

03 

-0.100G0E- 

-01 

w  • 

i     3  1  5  E 

0  3 

-0. 10000E- 

-01 

■  '  * 

1031 h ^ 

03 

C.22936E 

01 

V  • 

L  03 1 3  E 

03 

0.41 7  0 IE 

01 

V.    • 

1 0  3 1 2  E 

C  3 

0 .  c  5  E  6  r  E 

01 

o . ' 

■  }  2  r 

0  3 

3.89477E 

0  I 

'>.    • 

10312c 

03 

0. 10727E 

"2 

^.  • 

L  o  3  1  3  E 

03 

0  .  1  3  5  2  5  E 

0  ? 

"  3  i  ?  c 

03 

0.1 50 3  IE 

02 

_    • 

i .:  3 1 4- c 

0  5 

0.  177  83E 

(  2 

c  . 

L031  . 

0  3 

0. 1938  3E 

0  2 

0. 

L  I316E 

03 

0. 21607E 

-  2 

u  • 

L 3317E 

0  3 

0.  17  5  34E 

02 

.  . 

10  3  17  E 

u3 

0. 16  3  34E 

01 

0. 

10  31 6  E 

03 

-0. 10  000  E- 

-01 

J  • 

.0  315E 

03 

-0.100C0E- 

-01 

'-  . 

103 14E 

03 

0. 12404E 

30 

0. 

10 3 131 

0  3 

0.  198 12E 

01 

k.   . 

1 :  3  1  3  E 

u3 

0.*5373£ 

01 

L  0 31 2E 

03 

0.63256E 

01 

'.'. 

103 12E 

03 

0.3c 594E 

L  1 

t    . 

Li  312  E 

0  3 

0.11485E 

02 

1 . 

L0312E 

0  3 

0.  1303  IE 

02 

Ll  3  13L 

0  3 

0, 15  844E 

02 

'«_  . 

1 0  3 1 4  E 

0  3 

0. 1742  9E 

0  2 

0. 

L0315E 

03 

0 .  1 9 7 S  t  E 

0  2 

<J . 

I  J  3  1 6  £ 

03 

0  . 1 5  8  9  8  E 

0  2 

p. 

1 0  3 1  c  E 

03 

0.6] 559E 

00 

V.'   • 

L0315E 

03 

-0.100  00E- 

-  0  1 

p.: 

3314c 

0  3 

0. 78  560c 

00 

<j  • 

103 13E 

0  3 

0. 38016E 

01 

u  * 

L  0  3 1 2 1 

0  3 

0. 55  553E 

0  1 

0  • 

1 0  3 1 2  E 

03 

0.81511c 

01 

0. 

L 03 12c 

0  3 

0. 10233E 

0  2 

0. 

L0312E 

03 

0 .  1 2  1 7  3  E 

0  2 

0  . 

103 13E 

0  3 

0 .  1 4  8  4  7  £ 

02 

/"t 

.031 3E 

0  3 

0. 16352E 

02 

0. 

1 C  3  1  <»  E 

0  3 

0.19126 E 

0  2 

0. 

1<  315E 

0  3 

0.  17699E 

02 

C  . 

L0316E 

03 

0.58C 3  7 E 

(il 

c . 

L  0  3  1 :  :: 

0  3 

-o. iooooe- 

-01 

o  . 

10  31 4 E 

03 

0  .  6  3  V  9  3  b 

00 

c  • 

I  C  3  1  3  E 

03 

0.  3-f574E 

01 

r 

1C312E 

0  3 

0 . 51556 E 

01 
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HASAN    KJCAOGLU. SPEED    CONTROL    OF    A    DC  .  MCTGR  .  NO  .  59 


TIME 


VI  N 


VOUT 


IAR 


0.27522E  01 

0,275  3  1 

0  .  2  7  5  5  h-  t  .  1 

0  .  2  7  5  7 0  t  01 

0.27  5S5E  Gl 

0.276<   l  i 

0.276171  . 1 

0.276 33 E  ci 

0.27o4  i 

0.27665c  CI 

0.276606  CI 

D. 27  6' 6  6  CI 

0.277121  0  1 

0  .  2  7  7  2  6  E  01 

0.27744c  v.  1 

0 . 27760  £  C  1 

0.  2  11  li  E  CI 

0.27791c  Gl 

0.27  J u7  =  01 

0.27823c  v.  1 

0.273  39E  i  1 

0. 278 65 1  (  1 

u • 276  70  :  0  1 

0. 278  6c  ,1  i  1 

j  .  ^  ? i  •  ;  -  v  l 

0.279 16 i  01 

0  •  2  /  9  34 1  U 1 

0.2  79  5*  '  1 

0.27965 >  r  1 

U . 2  7 9  e 1 L  11 


0.10953E  03 

I . 113S4E  0  3 

0.10332c  C 3 

0.11500c  i'3 

0.10374E  03 

i  .  1 1 3  7  v.  t  0  3 

0 .  3 ** 3 5 3E  1--? 

0.256295  0  2 

0.10386E  03 

0.11 500 E  03 

0.1j319£  03 

0.11368!  l3 

0  .  1 0  9  4  8 E  03 

0.1103  3c  C3 

0.1 1345 E  (.3 

0  .  1 0  4 '-*■ 1  E  03 

0.1149  8c  03 

('.:^;;3l  02 

C.4«*276E  02 

0  .  <:  30  3  5 :-  G  2 

0 . 1 1  0  5  3 1  0  3 

0.11323!  03 

0. 1  J  6  >-,•  -.  Cj 3 

0.  11 49 3  c  03 

0.10  20  3c  03 

0.1142  IE  0  3 

0. 10663c  03 

0  .  1 1  1 0  * E  0  3 

0.11300E  03 


C  .  1 0  3 1 2  £  C  3 

0.10312b  03 

0.103126  03 

0  .  1 C  3 1  ^  l  0  3 

0.10313c  u3 

0.10  314E  2  3 

0.10315E  03 

C.1C316E  03 

C  .  10315c  03 

0.10 3 14E  0  3 

0.10  313E  0  3 

v.  103  13c  0  3 

Ci  .  1 0  3 1  2  E  0  3 
0  .  10  3 1 2 1 

0.10312c  u3 

0  .  1 0  3 1 2  E  C  3 

0.10313E  0  3 

G.  10314?  03 

0  . 1 1  315c  0 3 

0. 10316E  ):< 

0.  10315:  03    - 

i .10314c  03 

0.10313E  03 

U  .  1 0  3 1 3  E  0  3 

D  .  1 0  3 1 2  i  03 

0  .  1 0  3 1 2  J  0  3 

0.1031  2 E  0  3 

C. 103131  03 

0.103 13=  0  3 

0  .  10  3 1 4  i  0  3 


0.81 15 IE  ol 

0.97 685 E  01 

0.12  330E  02 

0.14354E  32 

0.162573  02 

0  .  1 8 7 95 E  0  2 

G.  193002  02 

0 .  1 G  S  0 3E  0  2 
■0.  100  00E-01 
•0. 10000E-01 

0.20634E  .-1 

0. 42040 E  01 

0.68  989E  01 

0.857036  01 

C.11476E  02 

0.13047E  0  2 

0.15679c  02 

0 . 1 7  3 1 3  E  02 

0.19 3 96E  02 

0.14003E  02 

0. 10000E-01 

0.29790E  00 

0. 20  3  3  3c  01 

0.48  705E  01 

0.6843  IE  01 

0.69  994E  01 

0 .  1 1  5  3  4 1  02 

0.1316 St  6  2 

C.15  936E  02 

0.1 7465 E  02 


0.27997 t    C 1 
0.23013c    01 


J.23021 


CI 


0.2 8046 E  01 

0.2 8060 6  01 

0.26076c  LI 

0.23092c  0 1 

0  .  2  6  i  (.  6  r  01 

0.25124E  CI 

C.  26  W  E  Gl 

0 . 2  3 1 5  5  E  01 

v  . .  2  6 1 7 1  £  .  ) 

0.231c 7 r  C 1 

0  .  2  6  2  0  3  £  0  1 

0 .  2  32 1 9  E  C  1 

0.2323*E  61 

0.2  82  50E  01 

O.ZoZii.i:  01 

0.262S2E  CI 

0.2629  6c  01 


C  .1  )546E  C  3 

0  *  5 1  4  6  p  E  ^  2 
C  .  C 
0  •  ^.; 

C .  10826E  03 

0.1113  5!:  03 

0  . 1 1 2  7  6  E  0  3 

0.  105971:  0  3 

0 .  1 1 4  6  5  E  0  3 

i- .  1  6  0 1  a  c  C  3 

0.11448E  6  3 

0.1u762E  C3 

0.11166E  0 3 

0.76390E  o2 

0.3  5  55  2c  0  2 
6 .  0 

0.10021c  u 3 

0.11459E  0  3 

C.10737E  0  3 

0  . 1 1 1 9  6  E  6  3 


0.103156 

03 

C. 20  084E 

C2 

0.103 16 E 

0  3 

0. 173622 

02 

0.163]-.: 

0  3 

0.35  648E 

01 

0.10  315E 

0  3 

-o.  looooe- 

-01 

C  .  1 0  3 1 4  E 

0  3 

0.15162c 

0  0 

G. 10313c 

0  3 

0.31660E 

01 

0. 1031 3E 

0  3 

0.48425c 

01 

0. 103 1 2E 

03 

0.7692  2E 

01 

0.  10312!. 

0  3 

0 .  9  6  5  2  c  E 

01 

0.10312E 

0  3  . 

0.11766E 

02 

0.10313c 

0  3 

0. 14161E 

02 

0.10313 E 

0  3 

0.15796E 

0  2 

U  .  1 0  3 1  > L 

03 

0. 1  S3  69  2 

0  2 

0.10315E 

0  3 

0.  16 0  66 - 

02 

G  .  1 0  3  1  c  6 

03 

0.77945E 

Jl 

0  . 1  0  3  1  5  E 

0  3 

-0. 100 00 E- 

-(.  1 

0.10316- 

0  3 

0.751  126 

0  0 

0.10313E 

0  3 

0.326 32£ 

01 

C.10313E 

03 

0.51 J27E 

01 

0  .  1 0  3 1  2  E 

0  3 

0. 60  1 346 

01 
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HASAN    KGCACGLU. SPEED    CONTROL    OF    A    DC. MOTOR. NO. 59 


TIME 


VI  N 


VOUT 


I  A* 


0  .  2  8  3 1 4  E 

0.28329!     0 
0 .  ?  c  ?  3    . 
0.28  561l 
a. 283  77      0 

C  .  t. :.  3  9  3  E 
0.28^0= 
2  V 
0  .26440 
0 . 2  84  56  £     0 

0.  2 34 72c    Ci 

0  •  2<   -   : 

0.26504 [ 

0. 23;  1-  . 

3.28  535!  0 
0.28551; 

0.28  5  67E  C 

0.  28  5  83 1  '" 

0 .2859^  v. 

0  .  2  8  6 1 4  l  C 

0. 2  86  3C£ 
0  .  2  8< 

0  •  28662  i  ■_ 
0.26678E  C 
0. 2 86 94 £  c 
0.2 £709 
.  2  3  i  2:  i 
0  .  2 87 4 1 E  C 
C.26757E  0 
0.26773c    0 

0.237 89 E    0 
0.288C4 
0.  28620 1 
0  .  2  8  8  3  6 1    0 

J .  2b 3  5 2 ..     ( 


2  0 


C    J       I. 


0.238 84 C  C 

0 .  2 3 3 99 1  C 

0  .  2  8  v  I  5  :.  0 

G.28931E  0 

0.28  947E  C 
0.28963 

0. 289796  C 

0 . 2  8 9  94  E  u 

C  .  2 90 1 0 £  C 
0.  29(  26 
0. 29042c 
0.29058E 


'.11223:  03 

G .10  6  9  3c  U 3 

C  .  1  1  k'Vt  0  3 

G. 996 luE  02 

0,ll46?fc  0  3 

C . 90  155 E  J  2 

>  •  3  «t  4  3  ^  c  02 

0.1 145 9£  0  3 

j .  1002 Vr  U3 

0.  11478c  ^  3 
0.10643b 

1 .  1 1 2  5 1 E  1 3 

0  . 1 1 1 5  3  E  I  3 

:   .  1  j  7  6  s  E  0  3 

C .  1 1 4  4  7  E  u  3 

0.1*  '3 

C  •  5  2  3  9  o  £  0  2 

i.O 

0.0 

0.11133E  0  3 

;.10c26£  13 

0  .  1  1  -t  3  <+  c:  0  3 
0.10146:!  03 

1  .11',  3  1  :.  0  3 

L.lv543!  03 

0.1130  I r  : 3 

0.111015  03 

0 .  1 0  8  (0 1  0  3 

0.&3811E  0 2 

0 .  5 1 1 3  7 1  01 
0.0 

0  .  1  0  4  9  2  E  (^  3 

0.11324E  13 

0.11067b  03 

0  .10911 L  u 3 

0.11  4 J  <+£  0  3 

0.102b 5 u  0 3 


0.11498i 


0  3 


0  .  1 0  4  3  9  E  0  3 

D  .  1 1 3  4 6  E  03 

r-.110il.  03 

0.10931$  03 

0.7083  9E  0  2 

0.76545c  01 

C  .  53133E  0  2 

C  .  1 0  3  3  4  E  L  3 


0.10312E  03 

0.10  312E  0  3 

0.10313E  0  3 

0 .  1 0  3 1  3  fc  0  3 

0.1031 4 E  0  3 

0  . 1 0  3 1  5  E  )  3 

0.10316E  03 
0  .  1 0  3  1 6  E 

0 . 1 0  3 1 5 1  03 

0  .  l 3  3  1  -'.  E  0  3 


0.96O71E  01 
0.12427E  02 
0 . 1 4 1 5 0  E  0  2 
0.16 3 97E  02 
0.16c  lie  02 
0.193 15 E  r:2 
0.J.2733E  02 
■0.  10000E-C1 
■0.100COE-01 
0.15  806E    01 


0 . 1 0  3  1 3  E 

0  3 

0.4C604E 

01 

0. 10312c 

03 

0.64366E 

01 

0.1031, 

0  3 

0.82 774E 

01 

0 .1031 

03 

0.11C 99E 

0  2 

0.1031 3E 

03 

0 .  12  6 4 4  E 

0  2 

0.103  1 

0  3 

0. 15443E 

0  2 

0. 10314c 

03 

0.170*  7  E 

0  2 

0.10315E 

0  3 

C. 19 37  IE 

02 

0  .  1 «   31  :  .- 

'J  3 

0  .  1  5  2 

j  2 

0  •  1 0  3 1 6  E 

03 

-  %  100 COE- 

-C] 

o .  i  o  315:; 

03 

-0.  10C00E- 

-  0  1 

: . . .  i  .  31 

0  3 

0.13863E 

01 

0 . 10  ■  i     E 

0  3 

0.  '>3  814£ 

01 

0 .  1  '„  3  I  2  t 

0  3 

0.ol543E 

0  1 

0.  10312; 

0  3 

0.  b6  5  -  6  £ 

01 

0.  1031  2E 

0  3 

0.  ]  0  Es6  K 

0  2 

0.10312E 

03 

0.127   •• 

0  2 

0  .  1  r>  3  1  3  i 

03 

0.  15  i 

j2 

U  .  10  3 1  \  ■■ 

0  3 

0.  16390E 

0  2 

0.10315c 

33 

0.  19647E 

0  2 

0.10316c 

03 

0.18035E 

0  2 

0.103 16 E 

0  3 

0.53505E 

01 

0 • 1 0 315 E 

03 

-C.  10000c- 

-01 

j  .  1  J  3 1 4  c 

0  3 

-0. 100 OOE- 

-  0  1 

0  .  1 0  3  1 3  E 

03 

0  ♦  2  5  7  5 2l 

0  1 

0.  103  13E 

0  3 

0.4  2  740E 

01 

0  .10  312c 

0  3 

0. 72420E 

01 

0.  103  1 2E 

03 

0 .89239 E 

01 

0.  10312E 

0  3 

0.  11505E 

)2 

0.  10  312E 

03 

0.13548E 

02 

0. 1  136  c 


U3 


0.29089E    01       0.10994E    03 


'). 10313-  03 

0.1  )314t  0  3 

C.10315E  03 

0.10316E  03 

0.10317c  03 

0.10313  £  3 3 

0.10317E  03 

0  •  1 0  3 1  (.  E  0  3 

0.10315E  03 

0.10 3 14E  0  3 


0.15416E  0  2 
0.180  HE  02 
G.19471E    02 

0.22197c  02 
0.20c32c  Kil 
O.^Ob95c  01 
-0. 10000E-01 
-0.10  M OE-01 
-0. 10000E-01 
-0.  100 D0E-01 
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HASAN    KOCAGGLU. SPEED    CONTROL    OF     A    DC.  MOTOR.  NO.  -9 


T  I  M  E 


VI  N 


VOUT 


IAR 


0.29105 E 

0] 

0 . 

,  1  J  9  3  *  E 

03 

0 

.  ln313E 

03 

0.3O2  02E 

01 

;  .29121 E 

01 

/  • 

,  ]  1 3  6  9 1 

0  3 

< 

.10313 E 

03 

0.473  90E 

01 

D  .  2  :<  1 

>1 

(  . 

1  3376c 

03 

■.  .■ 

.10312  :' 

03 

0.745 97 E 

01 

9.29  1: 

.'  1 

,115  ;  0  i 

03 

c 

.  1 0  3  1 2  E 

03 

0.94781E 

01 

0  .  2  9 1 1 

01 

►1032  ' 

03 

0 

.  1 0  3 1 2  E 

0  3 

0 .114 94E 

0  2 

0.291! 

CI 

0. 

,  1 1 3  6  5 1 

U2> 

1 

.103]  -  E 

0  3 

0. 140S0E 

02 

■J  .  2  r  2  r.  c 

CI 

0  ( 

,  1")  ,'37: 

03 

0 

.  1 0  3  1  3  E 

0  3 

0. 15  6  3  3E 

0  2 

0. 23^169 

CI 

u'  1 

.  3.  1j23E 

^' 

.  1  !-  3 1  4  E 

0  3 

0. 1S424E 

0  2 

0. 292  32= 

'  I 

72  jj2l 

0  2 

0 

.  ]   1 3  1  5  E 

0  3 

0. 17402E 

02 

0.292  s  3  - 

Jl 

0 

,  10212E 

0  2 

0 

.  1  ;313f- 

0  3 

0.61 5 92E 

01 

C .  2  9  2  6  \  - 

C  1 

' 

\J 

c 

.  1 0  3 1 5  E 

03 

-9.  10000E- 

-(•1 

0  .  2  9  2  7  9  t 

01 

1 

,102735 

0  3 

\-.' 

. 103 1h£ 

•J  3 

0. 146 97E 

01 

0.29295! 

.  i 

J  . 

,  1 1 4u  2 E 

0  3 

\. 

.19  3  13c 

0  3 

0.4 174  IE 

01 

0.29311: 

.  1 

r- 

, 10917c 

o  i 

u 

. 10312E 

03 

0.59010E 

01 

C  .  2  9  3  2  7  £ 

.'1 

3. 

►  1 1 U6  1 E 

0  3 

c 

.  i  j  3  i  2  9 

0  3 

3-.  SS44SE 

01 

0  .  2 93 4 3  E 

CI 

./  1 

1132  8c 

0  3 

c 

.  1  ■  3 1 2  E 

■j  3 

0 . 10 450 E 

H2 

0.29353E 

'  1 

••    1 

. 1046  4c 

03 

0 

.1031 39 

03 

0.131  V?t: 

3  2 

0.293  74i 

•-  i 

0  , 

,  1 1 4  9  6  E 

-.3 

f , 

.  1 0  3 1  3  E 

0  3 

).  14993E 

(J  2 

0. 29390  z 

C  1 

1  •.,  2 1  5  E 

03 

/  4 

. 103 14 £ 

'J  3 

0.17012E 

0  2 

j.2r   '■ 

Jl 

0< 

,  1 1 4 1  SE 

(    5 

■J 

.  i  c  3 1  ■:  - 

03 

0. 19425E 

0  7 

0.  294 22 E 

01 

i   , 

b  c  7  3  9  E 

02 

n 

.  1  0  3  1 6  E 

0  3 

0.20 2 13E 

02 

0.2;-  r  ;  .. 

01 

•  • 

292  96c 

02 

o 

.19  3  17  E 

0  3 

0.123   "   : 

3  2 

0  .  2  9  -?  5  3  E 

01 

0  . 

,  0 

w' 

.10316c 

0  3 

-0.100  0)  E- 

-01 

0.294 £9c 

0 1 

Q  , 

c 

.10315 E 

03 

-0.  10000c- 

-CI 

0. 294  3:>t 

01 

1  i  1 

. 11492c 

0  3 

c 

.  1 0  3 1 4  E 

0  3 

D. 96490E 

00 
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INPUT  C-^TA  RECORD 


ORDER  OF  EQUATIONS 
INITIAL  T  I  ■•  - 
FINAL  T  I  M  = 
STEP  S  I  Z.: 


0.2i?94E  01 
C.3500E  CI 
C  .  BOO  .-  E-04 


THE  NON-2     CONSTANTS i  C(  I)  ,  A  ! 

C (  1  )  =  0.1150c  )  3 

C(  2)    =  0.9  )00E-u2 

C(  3)  =  O.SuQc/c-Cl 

C (  4  )  =  0  . 1 3  8  U  E  01 

C(  5)  =  0.33     02 

C  (  6  )  =  C  .  1 0  0  J I     02 

C(  8)  =  O.ICO^E  02 

THE  NON-ZERO  INITIAL  CONDITIONS  ARE 

X  (  1  )  =  0  . 1 U  5  3  E  0  3 

X(  2)  =  0.7794c  00 


THu  COLLIN  HEADINGS  AND  THE  CCRRcSPONDING  VARIABLES  ARE 


TIKE 

V1N 

y  U! 
IAR 


X( 
X( 

X( 
X( 


)) 

5) 

1  ) 

2  ) 


THE  INDIVIDUAL  GRAPH  TITLES  AMD  THE  CORRESPONDING  VARIAriU'S 


OUT  ^UT  VS 
INPUT  VS 


]  '•"  . 


x(  ^)  vs.  x(  io) 

X(  5)  VS.  X(10) 
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HASAN    KOCAOGLU. SPccD    CONTROL    OF     A    DC . MCTDR. NO. 69 


TIME 


VI  N 


VOUT 


IAR 


0.259  39i  ( 

0.2595.5E  C 

C  •  2  5  9  7  1 E  0 

0.259  66  £  U 
G.260( 
0  •  2  6  0 1 8 e 
0.2603^ 

0.26C5f    .  l 
0  •  2  6  0  6  c  ; 

0.26081 1  L 

'.260  91 

0.26113  ■■:  J 

0.261292  0 

0  .  26  i  ;■  5  0 

0.26161!  C 

0.26176c  0 

0.26192c  0 
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